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Overview  of  ONR  sponsored  research  to  Dr.  B.  Greg  Mitchell,  and  long-range  goals 

A  thorough  understanding  of  marine  optics  is  a  mission  of  si^iificant  Naval  relevwce.  The 
ability  to  use  optical  measurements  (in  situ  and  remote^o  estimate  light  propagation  and 
phytoplankton  growth  rates  is  a  significant  objective.  This  objective  requires  a  detailed 
understanding  of  the  nature  and  variability  of  source  and  loss  terms,  and  an  ability  to  model  the 
relevant  rates.  Figure  1  summarizes  dominant  state  variables  and  transformations  occurring  in 
marine  optics,  and  the  contributions  by  the  principal  investigator  in  understanding  these 
processes.  Several  significant  problems  previously  have  not  been  resolved  with  respect  to 
measuring  or  modeling  the  processes  in  Figure  1.  The  long-range  objectives  of  the  principal 
investigator  is  to  continue  studies  of  the  various  components  and  rate  processes  shown  in  Figure 
1  in  an  effort  to  improve  our  understanding  of  marine  optical  properties  and  photo-physiology 
and  ecology  of  phytoplankton. 


Figure  1.  A  conceptual  model  of  optically  important  processes  in  the  ocean.  The  long-range  objectives 
of  our  research  include  a  detailed  understanding  of  the  processes  and  state  variables  specified  in  this 
figure.  The  state  variables  in  the  boxes  (S1-S6)  correspond  to  dominant  sources  of  optical  variability  in 
the  oceans.  The  principal  investigator  dedicated  bis  research  effort  to  understancfing  the  diverse 
aspects  of  this  system.  Methods  for  estimating  the  magnitude  and  characteristics  of  Si,  S2  and  S3  have 
included  analysis  of  particle  absorption  using  micropbotometry  (Iturriaga  et  al.,  1988);  macropbotometry 
(\UtcheU  and  Kiefer,  1984;  Mitchell  et  al.  1984;  MitcbeU  and  Kiefer,  1988a;  1988b;  Mitchell,  1990);  in 
situ  optical  profiling  resulting  in  discovery  of  regional  bio-optical  relationships  (Mitchell  and  Holm- 
Hansen,  1991a;  Mitchell,  19^);  a  novel  method  to  estimate  SI  fiom  moored  ra^ometers  (Stramski  et  al. 
1992);  and  aircraft  (Tanis  et  al.,  1990)  or  satellite  observations  of  biomass  (Mitchell  et  al.,  1991).  Optical 
based  models  of  phytoplankton  photosynthesis  for  laboratory  cultures  (Kiefer  and  Mitchell,  19^;  Swik 
and  Mitchell,  1991)  for  Antarctic  phytoplankton  (Mitchell  and  Holm-Hansen,  1991b)  and  for  Arctic 
phytopplankton  (Com  et  al.,  1992)  have  been  developed  to  model  {mcess  Kl.  Siegel  et  al.  (1989)  used 
an  in  situ  optical  approach  to  estimate  Kl  and  K2-»-IC5  for  populations  in  the  North  Pacific  Gyre.  Carder 
et  al.  (1991)  have  u^  optical  methods  to  partition  SI  fiom  S2  and  S3  for  remote  sensing.  Sedimentation 
rates  of  Antarctic  and  Arctic  phytoplankton  (K2)  have  been  studied  by  Mitchell  and  Holm-Hansen 
(1991b)  and  Wassman  et  ai.,  (1990). 


OBJECTIVES  FOR  GRANT  N00014-J-90-1052 


Tbe  objectives  of  the  woric  sponsored  under  ONR  grant  N00014-90<J-1052  was  to  apply  the 
absorption  methods  developed  during  FY89  under  grant  N00014>89-J-1071.  In  FY89,  the 
meth^ology  of  the  Quantitative  Filter  Technique  (QFT)  was  verified  and  published.  In  FY90, 
the  objectives  were  to  utilize  the  methods  in  the  study  of  specific  problems  in  marine  optics.  In 
particular,  we  sought  to  use  the  measurement  of  marine  particulate  absorption  to  enhance  our 
undeistaneing  of  the  bio-optical  properties  of  the  eastern  north  Pacific  observed  during  ODEX. 

The  main  focus  was  to  evaluate  the  bulk  inherent  optical  properties  of  Sl-i-S2-fS3  in  Figure  1, 
and  to  develop  indirect  optical  methods  to  partition  SI  from  S2+S3  using  the  spectral  reflectance 
of  the  ocean. 

TASKS  COMPLETED 

In  collaboration  with  Dr.  Kendall  Carder  and  Dr.  Raymond  Smith  (and  members  of  their 
laboratories)  a  detailed  study  was  undertaken  to  develop  a  model  of  the  relationship  between 
variability  in  pigment  specific  absorption  coefficients  of  marine  particulates  and  the  apparent 
optical  properties  of  relevance  for  remote  sensing  of  pigments.  The  modeling  approach  was 
validated  on  the  basis  of  the  observations.  For  many  remote  sensing  applications  concerned  with 
the  retrieval  of  phytoplankton  pigments  one  requires  a  separation  of  the  effects  of  detritus 
(S1+S2)  from  the  living  phytoplankton  absorption  (SI),  knowledge  of  the  q)ectral  shape  of  the 
chlorophyll  a  specific  absorption  coefficient  of  marine  particulates,  obtained  during  ODEX  using 
the  QIT  methc^  developed  by  the  PI,  allowed  improved  models  to  be  developed.  The  results  are  ■ 
now  published  (Carder,  et  al.,  1991) 

In  collaboration  with  Dr.  Ronald  Zaneveld,  a  weighting  function  has  been  developed,  which, 
when  applied  to  the  spectral  absorption  coefficients  determined  for  particulates,  provides 
spectrally  weighted  absorption  coefficients  which  are  directly  comparable  with  the  broad-band 
determination  of  beam  attenuation  coefficients  using  the  Sea  Tech,  Inc.  transmissometer.  The 
weighting  function  was  applied  to  samples  collected  during  ODEX.  The  fraction  of  particulate 
beam  attenuation  (at  665  nm)  attributable  to  particulate  absorption  is  variable  ranging  frdm  10- 
50%  depending  on  photoadaptation.  The  results  were  presented  at  the  NATO  Advanced  Study 
Institute  on  Particle  Optics.  A  manuscript  for  publication  is  in  preparation. 

In  collaboration  with  Dr.  Dariusz  Stramsld,  a  statistical  method  was  developed  to  derive  the 
spectral  diffuse  attenuation  coefficient  from  data  collected  from  a  single  spectroradiometer 
moored  in  the  ocean.  The  approach  required  determing  the  cross-correlation  coefficients 
between  the  surface  height,  varying  due  to  the  passage  of  waves,  and  the  variance  in  downward 
irradiance  at  the  single  moored  instrument.  Results  from  this  study  are  now  published  (Stramski 
etal.,  1992) 
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$  Estimation  of  downward  irradiance  attenuation  from  a  single 

moored  instrument 

Dariusz  Stramski,*  Charles  R.  BooiHt  and  B.  Greg  Mitchell$ 

(Received  14  Scpiemher  1V90;  in  revised  form  22  May  1991 ;  accepted  bjunc  1991 ) 

h  Abstract — A  method  for  estimating  the  diffuse  attenuation  coefficient  for  downward  irradiance. 

Afj(/).  from  time-series  data  of  irradiance  and  sea-surface  wave  elevation  is  presented.  Data 
acquisition  is  accomplished  with  a  fixed-depth,  single  moored  instrument  system.  Kj(/.)  is  obtained 
as  a  ratio  of  the  standard  deviation  of  natural  logarithms  of  irradiance  to  the  standard  deviation  of 
wave  elevation.  This  statistical  estimation  assumes  that  the  irradiance  fluctuations  arc  caused  only 
by  light  attenuation  over  a  fluctuating  vertical  path  length  associated  with  oscillations  of  surface 
wave  elevation.  The  limitations  and  the  possible  range  of  applicability  of  the  method  depend  on 
environmental  condiiions.  including  surface  wave  charactenstics.  water  clarity,  meteorological 
conditions  and  solar  elevation.  The  probability  density  and  power  spearal  analysis  are  useful  for 
interpreting  the  validity  of  statistical  estimates  of  Based  on  this  preliminary  concept,  some 
recommendations  for  an  improved  experimental  design  arc  proposed. 


INTRODUCTION 

The  diffuse  attenuation  coefficient  of  downward  irradiance.  Kj{X).  where  X  is  light 
wavelength,  is  among  the  most  important  physical  parameters  for  characterizing  the  upper 
ocean.  An  extensive  literature  includes  its  importance  in  primary  productivity  studies 
(Steeman  Nielsen.  1974;Kirk,  1983).  optical  classification  of  water  masses  (Jerlov.  1976; 
Smith  and  Baker,  1987b).  bio-optical  modeling  (Smith  and  Baker.  1978a;  Morel,  1988), 
remote  sensing  of  ocean  color  (Gordon  and  McClunev.  1975;  Gordon  and  Morel,  1983; 
Gordon  ei  al.,  1988)  and  considerations  of  heat  and  momentum  budgets  (Simpson  and 
Dickev.  1981;  Lewis  e/  al.,  1983;  Dickey,  1988).  Although  /Cd(A)  is  one  of  the  best- 
documented  optical  properties  of  the  ocean  (e.g.  Tyler  and  Smith.  197U;  Jerlov,  1976; 
Smith  and  Baker,  1978b;  Siegel  and  Dickey,  1987;  Morel,  1988).  our  understandingof  its 
variability  is  still  seriously  limited  because  the  data  typically  have  been  acquired  by- 
instruments  attached  to  a  cable  and  lowered  from  a  ship.  This  type  of  deployment  is  limited 
temporally  as  well  as  by  the  errors  in  estimate  of  near-surface  /Cu  due  to  wave  action 
(Weidemann  et  at.,  1990)  and  ship  shadow  (Helliwell  et  at.,  1990). 

Recent  development  of  optical  oceanographic  moorings  (Dickey  ei  al.,  1986;  Booth  ei 
al.,  1987)  promises  a  powerful  tool  for  long-term  environmental  monitoring  with  high 
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Fig.  I  Geometrical  conditions  for  statistical  determination  of  irradiance  attenuation  coefficient. 


lemporal  resolution.  A  variety  of  problems,  otherwise  inaccessible,  can  now  be  addressed 
using  data  collected  with  these  moorings.  The  question  addressed  in  this  paper  is 
principally  of  a  methodological  nature.  Our  objective  is  to  estimate  Akj(>l)from  time-series 
records  of  downward  irradiance.  £d(A),  and  surface  wave  elevation,  6z,  taken  simul¬ 
taneously  by  a  single  depth  moored  instrument.  We  were  motivated  by  observations  that 
variance  of  Ejf/)  is  similar  in  spectral  shape  to  (Booth  et  al.,  1987).  A  simple 

statistical  algorithm  is  proposed,  and  its  usefulness  and  limitations  are  evaluated.  Our 
approach  is  unique  in  that  it  is  based  upon  a  single  fixed-depth  measurement.  Traditional 
methods  require  either  two  instruments  deployed  at  different  depths  or  a  single  instrument 
operated  in  a  profiling  mode.  While  :n  traditional  measurements  wave-induced  fluctu¬ 
ations  of  £d(A)  represent  an  undesired  noise  that  must  be  eliminated  by  time-averaging, 
our  method  uses  these  fluctuations  to  determine  /fd(^)- 

DEVELOPMENT  OF  THE  METHOD 

/Cd(/)  is  defined  as  the  rate  at  which  the  natural  logarithm  of  downward  irradiance. 
In  £d(/.).  decays  with  depth  z  (e.g.  Preisendorfer,  1961): 

A.;d(z,  ;,)  =  -d(ln  £d(z.-l))/dz.  (1) 

where  z  is  positive  downward.  Our  procedure  for  estimating  from  a  fixed-depth 

mooring  (referred  to  hereafter  as  ‘statistical')  assumes  an  idealized  physical  situation  in 
which  fluctuations  in  arc  attributable  solely  to  light  attenuation  over  a  fluctuating  path 
length  associated  with  fluctuations  of  water  surface  displacement,  6z,  due  to  wave  motion. 
Geometry  of  measurement  can  be  imagined  in  which  this  situation  is  achieved  with  a 
reasonably  good  approximation  (Fig.  1).  This  involves  an  essential  simplification  regard¬ 
ing  the  motion  of  the  sea  surface. 

Although  ail  periods  arc  theoretically  present  in  the  wave  spectrum,  most  wave  energy  is 
usually  concentrated  in  a  narrow  frequency  band  (e.g.  Neumann  and  Pierson,  1966). 
These  most  energetic  components,  or  significant  waves,  are  of  particular  interest  to  the 
present  study.  Assume  that  as  the  significant  components  of  the  wave  spectrum  propagate, 
the  surface  moves  up  and  down  about  the  mean  level  and  is  nearly  flat  for  extensive 
horizontal  distances  and  hcncc  temporally.  Although  the  latter  assumption  is  never  strictly 
satisfled  under  natural  conditions,  our  approach  remains  tractable  if  a  portion  of  a  surface 
wave  crest  (or  alternatively  trough)  is  subtended  by  an  upward-looking  irradiance 
collector  within  the  angle  beyond  which  the  total  internal  reflection  takes  place,  i.c.  within 
Snell's  cone.  Thus  the  depth  z,  at  which  the  irradiance  collector  is  located,  cannot  be 
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arbitrarily  great  and  should  be  less  than  L/8.  where  L  is  the  length  of  a  significant  wave 
component.  This  estimate  is  based  on  a  tentative  assumption  that  the  linear  dimension  of 
the  surface  patch  within  Snell's  cone  does  not  exceed  L/4.  Typical  wavelengths  in  the 
ocean  range  from  about  20  to  200  m  (Defant,  196J).  Thus  the  permissible  depth  for 
measuring  irradiance  is  confined  to  the  uppermost  few  meters  if  relatively  short  waves 
prevail,  and  extends  to  about  25  m  when  significant  waves  become  very  long. 

Deployment  of  a  moored  instrument  at  such  small  depths  may  violate  our  fundamental 
assumption  that  variation  in  dz  is  the  dominant  source  for  producing  irradiance  fluctu¬ 
ations.  In  real  field  conditions,  the  geometric  structure  of  a  roughened  sea  surface  is  very 
complex  and  shorter  waves  are  superimposed  on  significant  carrier  waves.  Shorter  waves 
also  can  contribute  to  irradiance  fluctuations,  principally  at  relatively  small  depths  through 
processes  associated  with  fluctuations  in  surface  slope  and  curvature,  especially  under 
clear  sky.  At  depths  shallower  than  2Ka'  (a),  irradiance  fluctuations  can  be  dominated  by 
focusing  of  sunlight  beams  (Dera.  1970).  For  estimating  KaU)  the  moored  instrument 
should  not  be  deployed  within  this  layer  unless  the  focusing  effects  are  minimized  using  a 
special  design  for  the  irradiance  meter.  We  will  return  to  this  problem  below  when 
discussing  the  limitations  of  the  proposed  approach. 

Now,  we  assume  that  as  the  surface  undergoes  its  motions,  the  incident  irradiance  from 
sun  and  sky  does  not  change.  We  also  assume  that  is  invariant  both  in  time  and 
throughout  the  water  column  above  the  mooring.  This  appears  to  be  valid  over  the  time 
and  space  scales  in  question.  There  will  be  variations  in  underwater  irradiance  resulting 
solely  from  variations  in  water  column  height  above  the  collector.  We  do  not  investigate 
the  significance  of  the  mean  cosine  of  the  underwater  light  field.  If  the  surface  is  at  the 
mean  level,  the  associated  irradiance  is  ?.).  Shown  in  Fig.  1  are  two  other  positions  of 
the  surface  yielding  higher  (position  1 )  and  lower  (position  2)  irradiances  than  £d (z.  A).  By 
operational  definition,  /fo(A)  can  be  evaluated  at  the  surface  using  the  irradiance  readings 
at  the  mean  depth  z.  Taking  the  readings  corresponding  to  the  surface  position  1  at  level  dz 
and  mean  surface  level  0,  we  can  write; 


£d  =  £d(z  -  dz)  -  In  Ed(z)],  (2) 

where  A  has  been  omitted  for  the  sake  of  brevity.  The  identical  formula  can  be  written  for 
any  <5z  at  any  time  t  including  upward  surface  displacements  like  position  2  in  Fig.  1  (notice 
that  in  this  position,  as  opposed  to  position  l.dz  is  negative  and  In  £j(z  -  dz)  <  In  £d(z)). 
Taking  many  samples  over  time  it  is  possible  to  resolve  the  variations  of  both  In  Ej  about 
In  £d(z)  and  dz  about  its  mean  value,  which  isO.  Thus  ^vorugtwg  ovotime  T  anif  squaring 
each  side  of  equation  (2)^ancL  over  ■time.  T ; 


Kl  = 


ir 

tL 


dz(i)‘  di 


■'i 

T 


(In  £d(z  -  dz(/))  -  In  £u(z)]-  d/ 


(3) 


we  find; 


_a(ln£d) 


(4) 


where  a(ln  Ej)  is  the  standard  deviation  of  the  natural  logarithm  of  irradiance,  a{dz)  is  the 
standard  deviation  of  water  surface  displacement,  and  Ed(z  -  dz(/))  and  dz{i)  are  the 
instantaneous  values  for  downward  irradiance  and  water  surface  displacement  at  time  i, 
respectively.  An  alternate  statistical  derivation  for  estimating  from  water  column 
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height  and  irradiance  fluctuations,  suggested  by  Preisendorfer  (personal  communication, 
see  Booth  et  al.,  1987).  yielded  results  similar  to  those  presented  here. 

The  approximation  described  by  equation  (4)  makes  restrictive  assumptions  that  must 
be  considered.  The  ‘statistical'  estimates  of  will  be  inaccurate  to  the  degree  that 

sources  other  than  variations  in  dz  contribute  significantly  to  the  variance  of  A 

review  of  previous  work  on  irradiance  fluctuations  will  lead  us  to  criteria  which  must  be 
satisfied  for  reasonable  application  of  this  method. 

Fluctuations  in  underwater  irradiance  arc  associated  with  fluctuations  of  surface 
elevation,  surface  slope  and  surface  curvature  (Schenck.  1957;  Dera  and  Gordon.  1968; 
Snyder  and  Dera.  1970).  Typically,  multiple  surface  wave  processes  affect  instantaneous 
values  of  £j(/).  The  distinctions  among  these  processes  represent  an  extremely  difficult 
problem,  for  which  no  rigorous  solution  is  available.  However,  previous  simplified 
models,  as  well  as  experimental  data,  have  proven  useful  for  describing  general  features.  It 
is  appropriate  to  enumerate  briefly  features  that  are  of  particular  interest  to  this  study. 

First,  the  wave-induced  fluctuations  of  £d(A)  for  clear  skies  may  be  several  orders  of 
magnitude  greater  than  for  cloudy  skies  due  to  surface  curvature  focusing  of  direct  sunlight 
beams  (Schenck.  1957;  Snyder  and  Dera,  1970).  This  is  especially  true  at  depths  smaller 
than  £  J '  (Dera.  1970;  Dera  and  Olszewski,  1978;  Dera  and  Stramski,  1986;  Stramskj. 
1986a).  The  focusing  effects  vanish  under  completely  diffuse  irradiation  of  the  sea  surface 
(Dera  and  Olszewski,  1967;  Stramski.  1986b). 

Second,  focusing  is  dominant  within  the  high-frequency  portion  of  the  power  spectrum 
of  irradiance  fluctuations.  The  frequency  where  focusing  begins  to  dominate  is  inversely 
proportional  to  the  square  root  of  the  depth  (Snv  der  and  Dera,  1970;  Fraser  et  al. .  1980). 
The  effects  attributable  to  light  attenuation  associated  with  fluctuations  of  surface 
elevation  are  expected  to  be  dominant  at  low  frequencies  within  the  most  energetic  range 
of  the  wavc-elcvation  spectrum  (Snyder  and  Dera.  1970;  Nikolaev  era/.,  1972). 

Third,  as  the  short  waves  have  their  point  of  maximum  focus  nearer  the  surface  and  the 
longer  waves  focus  at  greater  depths  (Schenck.  1957;  Dera  and  Gordon.  1968;  Nikolaev 
cl  al. ,  1972),  the  high-frequency  components  of  irradiance  fluctuations  are  rapidly  damped 
with  depth,  causing  a  sharpening  of  the  pow  er  spectra  (Gordon  e/ of. ,  1971;  Prokopov  and 
Nikolaev.  1976). 

Fourth,  except  for  cloudy  skies,  the  contribution  of  focusing  to  total  variance  of  £^(7.) 
appears  to  be  negligible  only  at  depths  greater  than  2ACJ '  (if  measured  with  traditional  flat 
irradianee  collectors  having  a  surface  area  of  the  order  of  10  cm^  or  less).  In  this  lower, 
‘diffuse’  layer  of  the  photic  zone  the  focused  direct  light  is  of  little  importance  and  the 
magnitude  of  fluctuations  in  >s  comparable  to  that  in  the  upwelling  irradiance 

(Dera,  19711).  This  suggests  the  dominance  of  the  mechanism  associated  with  variations  in 
water  column  height. 

Finally,  according  to  the  single-ray  theory  (Snyder  and  Dera,  1970),  the  effects 
associated  with  focusing  and  variations  in  wave  elevation  are  the  only  important  mechan¬ 
isms  for  producing  irradiance  fluctuations  when  the  sun  is  high.  At  low  solar  elevation  the 
mechanisms  associated  with  fluctuations  in  surface  slope  (which  strictly  vanish  when  the 
sun  is  directly  overhead)  may  also  be  important,  including  fluctuations  of  path  length 
which  were  shown  to  be  linear  in  depth. 

These  observations  indicate  that  restrictions  in  the  usefulness  of  the  present  statistical 
method  will  oecur  under  clear  sky  conditions.  Wc  anticipate,  however,  that  the  method 
can  give  reasonable  estimates  of  /Cj(^)  if  the  depth  of  measurement,  z,  satisfies  the  criteria; 
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2/CJ'(A)S2<  L/8.  (5) 

If  this  holds  for  the  most  penetrating  light  wavelength  then  a  reasonable  accuracy  can  be 
expected  over  the  entire  visible  spectrum.  Equation  (5)  embodies  deployment  depth 
criteria  specifying  that  the  irradiance  meter  cannot  be  located  too  deep  because  too  large  a 
surface  area  would  be  seen  within  Snell's  cone,  nor  can  it  be  too  shallow  because  of 
focusing  effects. 

In  practice,  the  reliability  of  the  statistical  method  will  be  difficult  to  verify  by  equation 
(5).  The  length  of  significant  waves.  L.can  be  estimated  from  the  wave-elevation  spectrum 
and  the  surface-wave  dispersion  relation  but  is  unknown  a  priori  (in  fact,  it  is  to  be 
estimated).  Nevertheless,  equation  (5)  is  useful  because  it  indicates  that  the  statistical 
method  can  be  employed  when  >  16/Z,  (Fig.  2).  In  Fig.  2  the  line.  ATj  =  16/L. 

delimits  two  regimes  represented  by  surface  wave  (L)  and  water  clarity  (Ka)  parameters. 
Failure  of  the  method  (/^d(‘^)  <  16/Z,]  may  occur  for  clear  waters,  where  is  of  the  order 
of  10'“  m~‘.  and  for  slight  sea  states. 

However,  we  believe  that  equation  (5)  represents  conservative  criteria  for  the  accept¬ 
able  depth  of  measurement.  It  is  possible  that  the  approach  will  still  be  reliable  if  2  £  Z-/4. 
Furthermore,  it  seems  that  the  range  of  applicability  can  be  substantially  increased  if  the 
irradiance  meter  is  fitted  with  a  diffuse  collector  having  a  large  surface  area  of  the  order  of 
10'  or  even  10^  cm^.  Because  the  collector  is  essentially  the  spatial  integrator,  the  focusing 
effects  will  be  masked  over  that  large  surface  area.  Thus  the  criterion  z  a  2/f  can  be 
perhaps  safely  replaced  by  2  a  Af  if  the  collector  is  sufficiently  large.  Consequently,  a 
much  less  restrictive  overall  criterion  could  be  K^i  2  4/Z-;  the  corresponding  line  is  also 
shown  in  Fig.  2.  Another  point  is  that  under  cloudy  sky,  when  the  above  water  light  field  is 
completely  diffuse,  the  focusing  effects  and  other  unwanted  contributions  to  fluctuations 
of  £d(^-)  3re  vanishingly  small.  The  statistical  method  should  then  successfully  predict 
assuming  z  s  US,  or  less  conservatively  2  s  Z,/4. 

MEASUREMENTS  AND  COMPUTATIONS 

Data  were  acquired  during  the  November  1985-June  1986  deployment  of  a  mooring 
system  approximately  1000  m  seaward  of  the  Scripps  Pier  (32°52'N.  1  IT’IT'W).  A  detailed 
technical  description  of  this  project  is  given  elsew  here  (Booth  e/o/..  1987).  The  instrument 


Fig.  2.  Curves  representing  the  lower  boundaries  for  the  range  of  applicability  of  the  statistical 
method.  More  or  less  conservative  criteria  were  assumed  for  estimating  the  upper  and  the  lower 
curve,  respectively  (sec  text  for  details). 
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w;is  deployed  on  u  3-point  mooring  at  a  nominal  depth  of  8  m  below  the  surface  Actual 
depths  ranged  from  5  to  Mm,  accounted  for  by  the  tides  and  by  the  changes  in  flotation 
during  the  project.  The  time  variations  in  downwelling  irradiance.  £<)(/./).  v^ere 
measured  at  13  wavelengths  (A),  viz.  410.  441 . 488.  520.  .350.  560.  589.  633.  656.  671 . 683, 
694  and  710  nm.  by  a  moored  spectroradiometer  having  a  6.3  cm  diameter  cosine  collector 
(MER-1064.  Biosphcrical  Instruments.  Inc.).  Detection  over  narrow  spectral  bands 
(<10nm  half-bandwidth)  was  provided  by  a  combination  of  absorbing  and  interference 
filters.  The  time  variations  in  water  column  height.  r(/).  due  to  surface  gravity  waves  above 
the  moored  platform  were  measured  by  a  pressure  transducer.  No  corrections  to  the 
standard  hydrostatic  equation  were  made  when  converting  pressure  readings  to  depth, 
which  leads  to  negligible  underestimation  of  dz{i)  (Colli.ms.  1976). 

In  this  study  we  use  data  sampled  at  0.4  s  intervals  for  cither  6.8  or  13.7  min  for  a  total  of 
either  1026  or  2052  digital  data  points.  The  averaging  time.  T.  in  equation  (3)  equals  this 
record  length.  The  scan  time  required  for  sampling  all  discrete  wavelengths  and  depth  was 
10  ms.  allowing  a  ‘snapshot’  measurement  of  £j(A)  and  dz.  Given  the  record  length  and 
sampling  rate,  the  range  of  frequency  which  can  be  analysed  is  between  ca  10”*  Hz  and 
1.25  Hz.  This  range  encompasses  completely  the  energetic  range  of  wave-elevation 
spectra,  but  a  significant  magnitude  of  oscillatory  power  of  wave-induced  irradiance 
fluctuations  at  frequencies  >1.25  Hz  may  in  reality  occur  (Gordon  ei  al..  1971 ;  Dera  and 
Olszewski.  1978). 

In  addition  to  the  moored  instrument,  a  similar  vertical  profiling  system  was  used 
periodically  at  a  location  within  100m  of  the  moored  instrument.  The  /(d(A)  values 
referred  to  hereafter  as  'reference'  were  calculated  from  the  vertical  profiler  according  to 
equation  ( 1 )  over  the  depth  interval  from  2  to  8  m.  This  calculation  is  similar  to  that  used  by 
Smith  and  Baker  (1984)  and  Siecel  and  Dickey  ( 1987).  The  vertical  profiling  system  also 
included  a  transmissometer  for  measuring  the  beam  attenuation  coefficient  at  665  nm. 
w  hich  was  used  to  monitor  vertical  structure  of  the  optical  properties  in  the  near-surface 
layer. 

We  have  carried  out  computations  of  the  ‘statistical’  K.,  from  equation  (4).  using  data  on 
temporal  fluctuations  of  irradiance  fuft.  A)  and  water  column  height  z(t)  as  estimated  by 
the  moored  instrument.  The  system  deployment  and  data  sampling  were  not  specifically 
designed  to  test  the  proposed  method,  so  a  limited  amount  of  data  is  now  available.  A  total 
of  1 1  time-series  records  taken  on  three  dates  were  chosen  as  adequate  because  high 
sampling  frequency  was  possible,  the  platform  was  less  than  .5®  from  level,  and  ncar- 
simuitancous  skiff-based  vertical  profiling  occurred.  Of  these,  nine  were  taken  on  the 
same,  almost  cloudless  day.  The  changes  in  platform  orientation  caused  by  surface  wave 
motion  were  typically  only  within  ±3°  from  the  mean  angle,  and  have  a  negligible  effect  on 
irradiance  measurements  with  a  cosine  collector. 

The  raw  irradiance  values  £o(/)  were  log-transformed  and  tested  for  the  presence  of 
slowly  varying  trends  that  can  be  due  to  changing  atmospheric  lighting  conditions  (i.c. 
changes  in  solar  elevation,  passage  of  clouds).  The  run  test  was  applied  to  this  problem 
(Bendat  and  Piersoi..  1966).  If  no  underlying  trend  was  present,  the  log-transformed 
data.  {In  £j((f)),  were  transformed  to  a  new  time  history  record.  In  having  a  zero 

mean  value: 


ln£j(/)  =  {ln£,(0} 


'ln£j(0dt. 


(6) 
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The  reason  for  representing  the  original  data  by  {In  £d(/)}  instead  of  in  £.,(0  is  to  have  the 
In  £d(0  notation  indicate  a  zero  sample  mean  value.  When  the  mean  of  a  time  scries 
(In  £d(f)}  was  found  to  change  over  time,  then  a  transformation  was  needed  to  induce  a 
stationary  series.  This  was  done  by  subtracting  off  a  linear  trend  which  was  fitted  to  the 
data  by  a  least-squares  method.  The  term  comprising  an  integral  in  equation  (6)  was  thus 
replaced  by  In  £d(f)  as  predicted  from  the  regression  analysis.  A  series  transformed  this 
way  had  a  mean  value  of  virtually  zero.  In  cases  where  a  single  linear  function  did  not 
provide  adequate  representation  of  a  trend  over  the  entire  record  length,  a  time  series  was 
broken  up  into  parts  and  the  procedure  was  applied  to  each  part  separately. 

The  mean  was  also  removed  from  each  record  of  z{t)  to  obtain  a  time  series  of  wave 
elevation: 

<3z(f)  =  V^'(f)d/-z(/).  (7) 

T  jit 

No  significant  non-stationarity  was  identified  in  time  scries  z{i)  on  the  temporal  scales 
considered. 

Except  for  standard  deviations  needed  in  equation  (4).  we  have  estimated  probability 
density  functions,  power  spectral  density  functions  for  single  time-series  records  as  well  as 
cross-spectra,  and  hence  the  phase  and  coherence  functions  between  irradiance  fluctu¬ 
ations  and  wave  oscillations.  The  power  spectra  calculations  were  made  using  a  conven¬ 
tional  Blackman-Tukey  method  with  the  Parzen  lag  window  weighting  function.  The 
standard  error  is  22%  and  the  resolution  bandwidth  was  0.02.^  or  0.05  Hz  in  these 
calculations.  Full  details  of  the  computational  procedures  are  given  in  Bendat  and  Piersol 
(1966). 


RESULTS  AND  DISCUSSION 

Figure  3  shows  the  time  scries  of  the  downwcllihg  irradiance.  £j(/).  and  platform  depth. 
z(t),  recorded  on  21  March  1986  at  noon  (sun  elevation  56°.  clear  sky).  Maxima  in  E^U) 
coincide  w'ith  minima  in  z(i)  w'hich  suggests  that  we  arc  seeing  the  effects  of  the  path  of  light 
from  the  surface  to  the  instrument,  changing  with  the  surface  wave  motion. 
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Fig.  3.  Examples  of  lime  series  of  sea-surface  elcvaiion  and  spcciral  downward  irradiance  ai 
550  nm  measured  concurrent!)  hy  a  single  moored  instrument  under  clear  sky  conditions. 


8 


D  SiKAMSkI  (7  III. 


Figure  4  cxiimincs  the  probability  density  function  of  irradiance  data  presented  in  Fig.  3 
The  frequency  histogram  of  standardized  In  Ej(/)  with  a  mean  of  zero  and  a  standard 
deviation  of  unity  is  shown  along  with  a  fitted  normal  density  function.  The  x~  goodness- 
of-fit  test  indicated  that  the  observed  histogram  is  statistically  equivalent  to  a  normal 
density  function  at  a  level  of  significance  P  =  0.01 .  The  wave  data  recorded  by  pressure- 
sensing  devices  located  at  depths  of  interest  are  expected  to  represent  closely  a  Gaussian 
process  (Neumann  and  Pierson.  l%b).  Consequently.  In  £j(/)  is  expected  to  follow  a 
Gaussian  distribution  if  the  fluctuations  of  irradiance  arc  principally  caused  by  fluctuations 
of  d:.  The  Gaussian  oscillations  of  dz  will  simply  undergo  a  linear  transformation 
according  to  equation  (2).  so  that  the  output.  In  will  still  be  a  Gaussian  process.  On 
the  other  hand,  irradiance  fluctuations  dominated  by  focusing  effects  at  small  depths  are 
asymmetric  relative  to  the  mean,  due  to  the  non-linear  character  of  the  refraction  at  the 
air-water  interface  (DERAand  Stramski.  1986).  Whether  or  not  the  distribution  of  In  E^U) 
deviates  from  normality  is  important;  we  expect  the  statistical  method  to  be  invalid  if  the 
distribution  is  not  normal.  However,  the  question  whether  normality  of  ln£a(/)  is  a 
sufficient  criterion  to  accept  our  method  requires  a  better  experimental  design. 

The  power  sper'ra  of  In  £^(0  and  6z{t)  are  compared  in  Fig.  These  functions" were 
calculated  from  the  same  time-series  data  as  show  n  in  Fig.  3.  In  Fig.  5a  the  dominant  peak 
of  the  irradiance  power  spectrum  coincides  roughly  with  the  most  energetic  range  of  the 
wave  spectrum.  The  increase  in  the  cumulative  power  spectrum  of  irradiance  parallels 
closely  that  corresponding  to  surface  waves  (Fig.  5b).  These  features  again  suggest  that  the 
major  contribution  to  irradiance  fluctuations  is  due  to  fluctuations  of  surface  elevation. 

This  conclusion  can  also  be  supported  by  showing  that  focusing  of  light  cannot  account 
for  the  observed  maxima  in  the  irradiance  power  spectrum.  While  the  maxima  of  both 
spectra  cover  the  frequency  band  between  0.075  and  0.125  Hz.  the  dominant  frequency  of 
irradiance  fluctuations  due  to  focusing  at  a  platform  depth  (6.6  m)  should  be  above  0.7  Hz 
according  to  previous  observations  (Fraser  et  a!..  1980).  In  fact,  we  observe  very  little 
contribution  to  the  variance  of  irradiance  at  frequencies  higher  than  0.7  Hz.  In  addition, 
using  a  formula  which  describes  the  focal  depth  of  a  sinusoidal  wave  as  a  function  of  wave 
steepness  (the  ratio  of  wavelength  to  wave  height)  and  sun  elevation  (Khulapov  and 
Nikolaev.  1977).  we  estimated  that  the  significant  wave  would  give  rise  to  the  focal  depth 


STANOAROIZCO  In  e  4  I 

Kip.  4.  Kri-qucncy  hisitipram  of  ihc  niiliir;il  loparilhms  of  imidiuncc  dala  shown  in  Fig.  and  ihc 
hilcd  normal  prohahilily  dcnsiiy  funclion. 


MARKED  PROOF 
PLEASE  RETURN 


Estimalion  of  irradiancc  attenuation 


9 


0.00  0.29  0.90  0.79  t.OO  1.29 

fRCQUENCY  (H*) 


Fig.  5.  Estimates  of  power  spectral  density  (a)  and  cumulative  power  spectrum  (b)  for  wave 
elevation  and  natural  logarithms  of  irradiancc  data  shown  in  Fig  For  comparison,  the  power 
spectra  were  normalized  b\  the  total  variance.  The  cumulative  spectrum  at  any  frequency  / 
represents  the  fraction  of  total  variance  contained  between  0  and /. 


much  greater  than  the  platform  depth.  In  this  calculation,  we  assumed  that  the  significant 
waves  are  characterized  by  lengths  above  100  m  which  result  from  the  surface-wave 
dispersion  relation  =  Z-rg/L.  where  o/mj*  is  the  angular  frequency  at  the  maximum 
wave  spectrum  and  g  the  gravitational  acceleration).  The  height  of  these  major  wave 
components  is  approximately  4a(d:)  (e.g.  Massel.  1989).  which  yields  0.52  m  in  the 
present  case. 

In  Fig.  6a  the  coherence  values  near  unity  indicate  that  dz(i)  and  In  £d(0  ^re  highly 
correlated  at  low  frequencies  around  0.1  Hz.  In  a  first  approximation,  this  can  be 
attributed  to  a  physical  linear  system  with  a  constant  parameter,  relating  these  two 
variables  as  input  dz(i)  and  output  In  £u(f)  (e.g.  Bendat  and  Piersol,  1966).  Such  a  linear 
system  is.  in  fact,  the  es.scnce  of  the  derivation  of  the  present  method  for  estimating  Xd(<l) 
as  shown  in  equation  (2).  In  contrast  to  the  low-frequency  band,  the  variables  62(1)  and 
In  £d(i)  arc  almost  completely  unrelated  at  higher  frequencies  which  is  indicated  by  very 
low  coherence  values. 

A  function  shown  in  Fig.  6b  describes  the  phase  shift  between  62(1)  and  In  £d(0  at  any 
frequency.  These  variables  arc  nearly  180“  out  of  phase  at  low  frequencies  which  is 
consistent  with  the  shift  between  the  peaks  in  the  time  series  2(f)  and  £d(0  in  Fig.  3.  It  is 
natural  to  expect  that  such  a  shift  is  associated  with  light  fluctuations  induced  by  variations 
in  water  column  height.  The  observed  oscillations  of  phase  at  higher  frequencies  arc  likely 
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Fig.  6.  Estimates  of  coherence  (a)  and  phase  (b)  function  between  wave  elevation  and  natural 
logarithms  of  irradiance  data  shown  in  Fig.  3. 


Fig.  7.  Example  of  vcr>'  good  agreement  between  the  statistical  and  reference  estimates  of  the 
spectral  attenuation  coefficient  for  downward  iuadiancc. 


meaningless  because  the  magnitude  of  variance  of  both  processes  in  this  band  is  vanish¬ 
ingly  small. 

The  statistical  functions  discussed  above  bear  important  relationships  to  the  mechan¬ 
isms  of  irradiance  fluctuations.  The  benefits  of  these  functions  are  apparent  for  judging 
whether  or  not  the  statistical  method  for  estimating  Kj(A)  is  applicable  to  the  case  in 
question.  In  the  particular  example  discussed  above,  we  find  an  excellent  agreement 
between  the  statistical  and  reference  Ka(^)  (Fig.  7).  The  criteria  for  the  applicability  of  the 
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statistical  method  appears  to  be  nearly  satisfied;  the  platform  depth  was  somewhere 
betw'een  A'j‘  and  L/8,  and  the  statistical  properties  in  Figs  4-6  clearly  demonstrate  the 
dominance  of  fluctuations  of  path  length  associated  with  surface  elevation. 

Figure  8  illustrates  how  a  slowly  varying  trend  can  affect  the  present  method  for 
estimating  K^-  The  measurement  was  made  just  before  sunset  on  the  same  day  as  the 
previous  example.  The  trend,  apparently  caused  by  changes  in  sun  angle  under  a  clear  sky, 
is  evident  in  the  irradiance  time  series  (Fig.  8a).  As  shown  in  Fig.  8b.  the  statistical  method 
gives  erroneous  results  when  a  trend  is  present.  The  removal  of  the  trend  proves  effective 
in  improving  the  statistical  estimates  of  A'j(/).  The  detrended  data  provide  a  Aj(/) 
spectrum  which  bears  a  good  resemblance  in  shape  to  the  reference  spectrum  but  the 
magnitude  is  underestimated. 

A  summary  of  the  prediction  of  from  nine  measurements  at  13  separate 

wavelengths  made  on  21  March  1986.  is  included  in  Fig.  9.  These  data  correspond  to 
coastal  water  type  2  or  3  which  is  characterized  by  minimal  Aj  at  550  nm  of  about  0.2  m“' 
(Jerlov,  1976).  The  agreement  between  the  statistical  and  reference  values  for  is 
reasonable.  A  linear  regression  fit  provided  a  slope  of  0.97  and  an  intercept  of  -0.038  that 
are  not  significantly  different  from  1  and  0.  respectively  (0.01  <  P  <  0.05).  The  statistical 
estimate  of  A^.  on  average,  differs  by  17%  from  the  reference  value. 


TIME  (  min  ) 


Fig.  8.  (y)  Example  of  wavc-induccd  irradiancc  fluciuatinns  (at  .VSdnm)  superimposed  upon  a 
slowly  varying  trend  due  to  changing  sun  angle,  (b)  The  statistical  estimates  of  the  spectral 
attenuation  coefficient  based  upon  the  non-staiionary  and  detrended  time-series  data  compared 

with  the  reference  spectrum. 


Fiji  9.  The  statistical  estimates  of  the  irradiance  attenuation  eoeffieieni  plotted  vs  the  reference 
estimates.  Included  are  I  lf<  data  points  from  nine  separate  eompahsons  at  13  light  uavelengths. 
The  straight  line  represents  perfect  agreement  between  the  two  estimates. 


A  number  of  factors  that  are  not  related  to  the  basic  assumptions  of  the  proposed 
method,  may  contribute  to  the  observed  discrepancy.  First,  there  is  inaccuracy  associated 
with  measurements  with  both  the  moored  instruments  and  the  profiler.  Second,  the  depth 
measured  by  means  of  a  pressure  transducer  is  not  precisely  related  to  the  vertical  water 
path  below  a  wavy  surface.  Third,  changes  in  were  observed  with  the  vertical  profiler 
during  the  day.  Because  the  data  acquisition  by  the  moored  and  profiling  systems  was  not 
simultaneous  (displacement  between  the  compared  measurements  never  exceeded  1  h) 
and  profiles  were  made  50-100  m  from  the  mooring,  some  discrepancies  between  the  two 
estimates  are  expected.  Fourth,  some  data  were  acquired  at  low  solar  elevation.  This  may 
reduce  the  accuracy  of  the  statistical  estimation  as  suggested  earlier.  Finally,  the  water 
column  above  the  platform  w  as  inhomogeneous  as  observed  with  a  beam  transmissometcr 
(the  transmittance  per  meter  at  665  nm  increased  by  several  percent  between  the  surface 
and  the  platform  depth).  This  may  be  a  possible  source  of  the  observed  discrepancies 
because  the  reference  estimates  of  A'j(A)  are  based  upon  vertical  measurements  of 
irradiance  between  2  and  8  m  depth. 

Below  is  a  brief  discussion  of  two  other  examples  that  may  potentially  pose  problems  for 
application  of  the  proposed  method.  Figure  10  examines  the  overcast  sky  conditions.  In 
(a),  the  two-fold  change  in  irradiance  due  to  changing  cloud  coverage  on  the  scale  of 
minutes  is  shown.  Superimposed  upon  these  slow  variations  arc  short-term  fluctuations 
caused  by  surface  w’avcs.  ITtc  statistical  method  for  estimating  as  applied  to  the 

original  time  scries  of  irradiance.  is  obviously  invalid  (Fig.  10b).  The  correction  for  linear 
trend  w'as  applied  to  four  segments  of  the  In  £a(t)  series.  With  the  trend  removed,  this 
method  still  has  only  limited  success  because  of  the  poor  quantitative  agreement  with  the 
reference  estimation.  In  the  case  when  the  light  field  is  supposed  to  be  eompletely  diffuse, 
this  disagreement  appears  rather  unexpectedly.  We  have  no  evidence  for  any  specific 
explanation,  but  some  of  the  reasons  mentioned  earlier  may  be  responsible  for  the 
observed  discrepancy.  Wc  note  that  the  measurements  with  the  moored  instruments  were 
delayed  by  about  3  h  relative  to  the  vertical  profile. 

Figure  1 1  plots  time  series  of  wave  elevation  and  irradiance  measured  on  a  sunny  day  for 
relatively  clear  water  with  a  minimal  A'jof  about  0.1  m“'  at  488  nm  (this  falls  between 
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Fig.  10.  <a)  An  example  of  wave-induced  irradiance  fluctuations  (at  550  nm  l  superimposed  upon  a 
slowly  varying  pattern  due  to  cloud  coverage  changes:  13  February  i9S6.  sky  overcast,  sun 
elevation  .35°.  significant  wave  height  0.8  m  (h)  Comparison  of  the  statistical  estimates  of  A^(/.) 
before  and  after  trend  removal  with  the  reference  spectrum 


Fig.  1 1 .  Examples  of  time  scries  of  wave  elevation  and  downwiard  irradiance  at  4K8  nm  measured 
in  relatively  clear  water  under  clear  sky  conditions:  H  January  1986.  sun  elevation  28°,  significant 

wave  height  U.48  m. 
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Jerlov's  oceanic  water  types  II  and  III).  The  most  striking  difference  between  tltc  time 
scries  is  that  higher  frequencies  exist  in  the  irradiance  signal  as  illustrated  by  the  power 
spectra  (Fig.  12).  It  is  obvious  that  irradiance  fluctuations  at  higher  frequencies  (>b  3  Hz) 
cannot  be  associated  with  variations  in  surface  elevation.  Therefore,  a  simple  low-pass 
digital  filter  in  the  form  of  5-point  moving  average  (Fig.  13)  was  applied  to  the  irradiance 
time  series.  After  filtration,  the  higher  frequency  portion  of  the  irradiance  spectrum  down 
by  two  or  three  orders  of  magnitude  represents  negligible  contribution  to  the  total 
variance.  There  is  also  good  agreement  between  the  cumulative  spectra  of  waves  and 
irradiance  after  filtration  (Fig.  12b). 

The  statistical  estimates  of  A.'ij(A)  are  compared  to  the  reference  in  Fig.  14.  E\en 
though  the  undcsired  higher  frequency  fraction  of  variance  was  filtered  out,  the  statistical 
method  yields  a  considerable  overestimation.  This  proves  that  a  real  oceanic  situation  can 
result  in  failure  of  our  simplified  concept,  which  assumes  that  low-frequency  irradiance 
fluctuations  are  dominated  by  variations  in  water  column  height.  It  thus  appears  that  a 
significant  contribution  to  the  variance  of  irradiance  at  low  frequency  (ca  0.1  Hz)  ma>. 
under  certain  conditions,  be  attributable  to  fluctuations  of  surface  slope  or  curvature.  In 
the  present  case,  the  sun  elevation  was  low  (—28°);  as  mentioned  earlier,  this  can  result  in 
larger  effects  due  to  wave  slope  Furthermore,  whenever  wave  amplitudes  are  small,  as  in 
this  example,  the  relative  significance  of  changes  in  surface  elevation  to  the  variance  of 
irradiance  is  diminished. 


0.00  0.25  0.50  0.75  1.00 
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Fig.  12.  Eslimaics  of  power  spectral  density  (a)  and  cumulative  power  spectrum  (h)  for  wave 
elevation  and  natural  fogarithms  of  irradiance  data  shown  in  Fig.  II.  TTie  dashed  lines  represent  the 
irradiance  data  after  low-pass  frequency  liltration. 
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Fig.  13.  The  frequency  response  function  of  the  filter  which  was  applied  to  time-series  irradiancc 
data  and  yielded  the  dashed  line  spectrum  in  Fig.  12.  This  function  is  the  ratio  of  two  power  spectral 
density  functions  represented  by  the  dashed  and  solid  lines  in  Fig.  12. 


CONCLUSIONS 

Operational  and  experimental  mooring  deployments  will  certainly  continue  in  the 
future,  and  ultimately  benefit  various  long-term  oceanographic  programs.  The  proposed 
statistical  method  for  estimating  has  advantages  because  it  may  often  be  preferred 
or  possible  to  deploy  only  a  single  moored  system  rather  than  an  instrument  array 
configured  for  measurements  at  multiple  depths,  especially  if  cost  is  a  consideration.  With 
the  available  data,  the  method  has  some  limited  success.  While  reasonably  good  results 
were  obtained  when  the  actual  spectral  values  of  Kj  exceeded  0.2  m“',  the  failure  was 
observed  in  somewhat  clearer  waters.  The  method  appears,  however,  to  successfully 
predict  the  shape  of  the  Kj{/.)  spectrum  regardless  of  conditions.  This  is  an  important 
feature  because  our  approach  seems  to  be  applicable  even  before  further  improvements 
arc  achieved,  provided  that  A',|  is  known  for  at  least  one  wavelength. 

Our  method,  by  integrating  in  the  surface  layer,  may  have  greater  applicability  in 
coastal  or  continental  shelf  waters  which  typically  arc  more  turbid.  Such  waters  have 
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Fig.  14.  The  stalisiic:il  csiimaics  of  obiainvd  fn>in  original  irradiancc  data  and  after 

low-pass  frequency  filiruiion  compared  wiih  fhc  reference  spectrum. 
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shallow  depths  for  remote  sensing  signals,  and  it  is  difficult  to  estimate  using  traditional 
profilers  in  this  near-surface  region. 

The  results  of  this  study  can  be  used  as  a  framework  for  further  investigations  of  the 
proposed  method.  More  and  better  data  arc  needed  to  improve  our  understanding  of 
possible  restrictions  that  apply  to  the  method.  To  further  test  the  method,  the  experiments 
should  he  designed  to  include  accurate  measurements  of  surface  waves  and  irradiance 
fluctuations  under  various,  carefully  controlled  environmental  conditions.  The  measure¬ 
ments  fvir  calm  seas  and  low  sun  elevation  may  be  meaningless  for  this  analysis  because  the 
fundamental  assumptions  of  our  approach  will  not  be  satisfied.  Further,  a  large-surface 
irradiance  collector  or  multiple  irradiance  collectors  are  highly  recommended  in  order  to 
minimize  focusing  effects  that  prevail  in  the  near-surface  layer  under  clear  skies.  With  such 
a  design  a  reliable  decision  as  to  w  hether  or  not  the  statistical  approach  is  applicable  can 
probably  be  based  on  criteria  combining  the  depth  of  measurement  and  characteristics  of 
significant  components  of  the  surface  wave  spectrum.  The  mooring  instrument  should 
sample  and  process  data  so  that  sufficient  information  for  such  a  decision  is  stored  in  the 
data  record. 

It  is  proposed  that  real-time  statistical  analyses^-wrio^iwg  pow'cr  spteiFal  BwliiiuiDl  be 
performed  and  the  results  of  these  analyses  stored,  rather  than  storage  of  lengthy  high- 
frequency  data  records.  Better  understanding  of  relationships  between  the  dominant 
source  mechanisms  for  irradiance  fluctuations  and  their  statistical  properties  is  essential 
for  optimization  of  data  sampling  and  storage  procedures.  Therefore,  future  research  into 
the  applicability  of  the  method  requires  initial  studies  at  a  higher  sampling  ra»e  than 
presently  used  to  better  resolve  the  actual  irradiance  fluctuations.  In  such  a  case  the 
filtration  of  the  undesired  high-frequency  portion  of  total  variance  can  be  realized  by 
means  of  digital  techniques.  We  expect  that  such  a  study  would  provide  the  optimal 
sampling/proccssing  protocol  for  operational  deployments. 
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Reflectance  Model  for  Quantifying  Chlorophyll  a  in  the  Presence 
of  Productivity  Degradation  Products 
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Marine  colored  dissolved  organic  matter  (CDOM)  lalso  gilvin  or  yellow  substance)  absorbs  light  at 
an  exponentially  decreasing  rate  as  a  function  of  wavelength.  From  -HO  nm  to  about  640  nm. 
particulate  phytoplankton  degradation  products  including  pheopigments.  detritus,  and  bactena  have 
absorption  curves  that  are  similar  in  shape  to  that  of  CDOM.  In  coastal  areas  and  areas  downstream 
from  upwelling  regions,  these  constituents  of  seawater  often  absorb  much  more  light  than  do  living 
phytoplankton.  leading  to  errors  in  satellite-denved  chlorophyll  estimates  as  high  as  1 33*^.  Proposed 
N.ASA  sensors  for  the  1990s  will  have  spectral  channels  as  low  as  412  nm.  permitting  the  development 
of  algorithms  that  can  separate  the  absorption  effects  of  CDOM  and  other  phytoplankton  degradation 
products  from  those  due  to  biologically  viable  pigments.  .A  reflectance  model  has  been  developed  to 
estimate  chlorophyll  a  concentrations  in  the  presence  of  CDOM.  pheopigments.  detntus,  and  bactena. 
Nomograms  and  lookup  tables  have  been  generated  to  descnbe  the  effects  of  different  mixtures  of 
chlorophyll  a  and  these  degradation  products  on  the  R(4I2)  :R|443)  and  R(443)  :R{565)  remote¬ 
sensing  reflectance  or  irradiance  reflectance  ratios.  These  are  used  to  simulate  the  accuracy  of 
potential  ocean  color  satellite  algorithms,  assuming  that  atmosphenc  effects  have  been  removed.  For 
the  California  Current  upwelling  and  offshore  regions,  with  chlorophyll  a  s  1.3  mg  m"\  the  average 
error  for  chlorophyll  j  retnevals  denved  from  irradiance  reflectance  data  for  degradation  product-nch 
areas  was  reduced  from  r6l%  to  =23*^  by  application  of  an  algonthm  using  two  reflectance  ratios 
(R(4I2)  :R(443)  and  R(443)  :R(565))  rather  than  the  commonly  used  algorithm  applying  a  single 


reflectance  ratio  (R(443)  :R(56S)). 


1.  Introduction 

The  ne.xt  generation  of  ocean-viewing,  visible/near- 
infrared  satellite  spectrometers  will  include  for  the  first  time 
a  short-wavelength  band  which  will  permit  improved  chlo¬ 
rophyll  a  concentration  determinations.  The  sea-viewing 
wide  ficld-of-view  sensor  (Sca-WiFS).  the  high-resolution 
imaging  spectrometer  (HIRIS).  and  the  moderate-resolution 
imaging  spectrometer  (MODIS)  planned  for  launch  in  the 
1990s  will  all  contain  a  wavelength  band  in  the  410-420  nm 
range.  On  Sea-WiFS.  this  band  and  a  band  at  490  nm  will 
complement  other  bands  (443.  520. 565.  and  665  nm)  that  are 
similar  to  those  of  the  coastal  zone  color  scanner  (CZCS).  It 
will  also  have  near  infrared  bands  (765  and  865  nm)  useful  for 
separating  water-leaving  radiance  from  radiance  scattered 
from  atmospheric  aerosols  for  coastal  or  other  waters  with 
high  in-water  backscattering. 

Colored,  dissolved  organic  matter  (CDOM)  (also  gilvin  or 
yellow  substance),  pheopigments.  detritus,  and  associated 
bactena  (collectively  referred  to  hereafter  as  "degradation 
products.”  or  DP)  absorb  more  strongly  in  the  4l2-nm  band 
than  at  longer  wavelengths  (see  Bricaud  et  al.,  1981.  Kiefer 
and  Soohoo,  1982;  Kirk.  1976.  1980.  1983;  Kishino  et  al.. 
1985;  Mitchell  and  Kiefer,  1988n;  Roesler  et  al..  1989). 
whereas  phytoplankton  absorb  more  strongly  at  443  nm  than 
at  412  nm  [Jeffrey,  1980;  Bricaud  et  al..  1983).  These 
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contrasting  absorption  characteristics  can  be  exploited  for 
the  purpose  of  quantifying  viable  phytoplankton  pigments 
(mainly  chlorophyll  a  but  including  some  accessory  pig¬ 
ments)  independently  from  DP.  They  also  provide  a  means 
to  remotely  assess  the  absorption  due  to  DP. 

CZCS  algorithms  for  estimating  chlorophyfl  a  (Chi  o)  plus 
pheophytin  a  (pheo  a)  concentrations  perfonn  «|uite  well  for 
regions  of  the  ocean  where  scattering  and  absorbing  compo¬ 
nents  of  seawater  covary  with  these  pigments  [Gordon  and 
Morel.  1983;  Cordon  et  al.,  1983).  A  number  of  empirical 
and  semianalytical  optical  models  have  been  developed  to 
simulate  the  behavior  of  the  underwater  light  field  for  such 
Morel  case  1  waters  [Morel  and  Prieur,  I9T7;  Smith  and 
Baker.  1978;  Baker  and  Smith.  1982;  Gordon  et  al..  1988; 
.Morel.  1988;  .Mitchell  and  Holm-Hansen.  1991;  Mitchell. 
1991).  Such  models  have  been  used  as  the  basis  for  classi¬ 
fying  water  types  and/or  for  developing  remote  sensing 
algorithms. 

The  accuracies  of  these  models,  however,  may  decrease 
when  environmental  conditions  depart  from  those  represen¬ 
tative  of  the  data  set  used  to  empirically  derive  the  covari¬ 
ance  relationships.  For  instance.  DP  incltiding  CDOM  are 
produced  when  grazing,  photolysis,  and  other  mechanisms 
degrade  the  viable  plant  matter  at  and  dowiistream  from 
phytoplankton  blooms.  The  DP-to-chloropbyll  ratio  will 
change  dramatically  for  a  parcel  of  upwelled  water  over  a 
relatively  short  time,  from  chlorophyll-rich  and  DP-poor  to 
DP-rich  and  chlorophyll-poor.  Solid  evidence  for  the  occur¬ 
rence  of  this  scenano  can  be  found  in  two  separate  studies. 
Peacock  et  al.  [1988]  found  that  absorptioa  attributed  to 
CDOM  at  440  nm  was  at  least  sixteenfoU  that  due  to 
phytoplankton  pigments  within  an  offsboro  jet  from  an 
upwelling  region,  whereas  pigments  were  the  dominant 
absorption  agents  at  the  upwelling  center  near  the  coast. 
Similarly.  Carder  et  al.  [  1989)  found  that  pacticufaite  absorp¬ 
tion  at  440  nm  decreased  thirteenfold,  while  CDOM  absorp- 
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tion  at  440  nm  increased  by  60^  in  10  days  as  they  tracked 
a  phytoplankton  bloom  from  the  Mississippi  River  plume  to 
Cape  San  Bias.  This  widely  varying  CDOM-to-chlorophyll 
ratio  has  a  profound  effect  on  upwelled  radiance  in  the  blue 
443-nm  band  of  the  CZCS  and  a  smaller  but  still  significant 
effect  in  the  green  520-nm  band.  The  correspondence  in 
absorption  at  443  nm  and  520  nm  between  CDOM  and 
chlorophyll  creates  erroneously  high  estimates  of  pigment 
concentration  in  those  models  which  rely  solely  upon  either 
of  these  spectral  bands  to  indicate  absorption  due  to  phyto¬ 
plankton. 

Separate  studies  performed  in  the  Southern  California 
Bight  point  to  another  way  to  view  the  problem.  Smtih  and 
SuAer  [I982I  developed  a  pigment  algonthm  for  tnat  site  that 
appears  to  require  a  larger  absorption  etfect  per  unit  of 
chlorophyll-like  pigment  than  do  most  other  case  I  algo¬ 
rithms  [see  Gordon  and  Morel.  19831.  However,  the  algo¬ 
nthm  parameters  for  the  Southern  California  Bight  study- 
region  changed  with  time  (compare  Smith  and  Baker  [19821 
to  Smith  and  Wilson  [19811),  requiring  contemporaneous 
ship  and  satellite  data  for  parameter  adjustment.  Mitchell 
and  Kiefer  [1988«1  found  significant  vanability  in  pigment- 
specific  particulate  absorption  in  this  region,  which  may 
account  for  some  of  the  difference  between  the  two  algo¬ 
rithms.  but  the  presence  of  CDOM  and  other  degradation 
products  can  also  cause  such  variation.  In  any  event,  with¬ 
out  e.xpensive  surface  verification,  waters  laden  with  DP 
cannot  be  unambiguously  identified  based  upon  satellite  data 
from  sensors  using  as  limited  a  suite  of  spectral  channels  as 
the  CZCS  employed. 

Another  problem  is  that  models  which  do  not  separate 
absorption  due  to  viable  pigments  from  absorption  due  to  DP 
have  less  utility  to  researchers  applying  physiological  pn- 
mary  production  models  that  are  driven  by  quanta  absorbed 
by  photosynthetic  pigments  [e.g.,  Kiefer  and  Mitchell.  1983; 
Platt.  1986;  Smith  et  a!..  1989]  or  to  researchers  interested  in 
the  DP  rather  than  the  pigments. 

.A  previously  developed  semianalytical  reflectance  model 
is  modified  here  to  address  the  concerns  listed  above.  The 
model  is  used  to  develop  an  algonthm  that  utilizes  a  412-nm 
spectral  channel  in  addition  to  the  443-nm  and  565-nm 
channels  e.xpected  on  Sea-WiFS  in  order  to  estimate  from 
irradiance  reflectance  data  Chi  a  concentration  and  the 
absorption  effects  due  to  DP.  In  situ  reflectance  data  from 
the  California  Current  region  is  used  to  test  the  algorithm. 

2.  .Model  and  .Algorithm  Development 

Morel  and  Gordon  [1980)  describe  three  approaches  to 
interpret  ocean  color  data  in  terms  of  the  in  situ  optical 
constituents;  empirical,  semiempirical.  and  analytical.  In  the 
analytical  approach.  (1)  radiative  transfer  theory  provides  a 
relationship  between  upwelling  irradiance  and  the  in  situ 
constituents,  then  (2)  constituent  concentrations  are  denved 
from  irradiance  values  measured  at  several  wavelengths  by 
inversion  of  the  resultant  system  of  equations.  The  algorithm 
presented  here  uses  this  approach,  and  the  term  "semiana- 
lytical"  is  invoked  because  pieces  of  the  radiative  model  are 
expressed  by  empirical  relationships.  This  analytical  deriva¬ 
tion  of  the  two  quantities  in  question.  Chi  a  and  degradation 
products,  avoids  the  problem  of  statistical  nonvalidity,  due 
to  the  spatial  and  spectral  correlation  between  the  two.  that 
Prieur  and  Sathyendrunath  [1981)  encountered  in  their  at¬ 


tempts  to  decouple  the  optical  effects  of  those  two  quanti¬ 
ties.  In  the  following  sections,  the  radiative  transfer  model 
will  be  developed,  the  individual  terms  in  the  resultant 
equations  will  be  expressed  in  terms  of  Chi  a  and  degrada¬ 
tion  product  concentrations,  and  the  algorithm  for  determin¬ 
ing  the  constituent  concentrations  will  be  explained. 

2. 1 .  Reflectance  Mode! 

The  concept  of  using  spectral  ratios  of  water-leaving 
radiance  to  determine  phytoplankton  pigment  concentra¬ 
tions  has  enjoyed  widespread  currency  ever  since  the  idea 
first  surfaced  in  the  work  by  Clarke  et  al.  [I970|.  The  most 
commonly  applied  algonthm  for  use  with  CZCS  data  cor¬ 
rected  for  atmospheric  effects  falls  into  the  empincal  cate¬ 
gory  as  descnbed  by  .Morel  and  Gordon  [1980]  and  is 
descnbed  by 

[C]  =  A(V«  (I) 

where  [C]  =  [Chi  a]  -*■  [pheo  n)  (square  brackets  here 
indicate  concentration  in  mg  m"').  .A  and  B  are  empirical 
constants,  and  r,j  is  the  spectral  ratio  of  water-leaving 
radiance,  defined  as 

r,j  =  (2) 

[Gordon  and  Morel.  1983).  Symbols  not  defined  explicitly  in 
the  text  can  be  found  in  the  notation  section. 

Equation  (2)  derives  from  earlier  research  where  rep¬ 
resented  R{\,)lRi\j).  the  irradiance  reflectance  ratio.  Irra¬ 
diance  reflectance  is  defined  as 

/?' V I  =  £,(A,  0-)/Ej(A.  0-1  (3) 

[Morel  and  Prieur.  1977),  and  a  similar  quantity,  remote¬ 
sensing  reflectance,  is  defined  as 

R, , i \)  =  irL^lki/E jiA.  (4) 

[see  Austin.  1974;  Carder  and  Steward.  19851,  where  £„(A. 
0-)  is  upwelling  irradiance  just  below  the  sea  surface  and 
£j(A.  0-)  and  £j(A.  0-v)  are  downwelling  irradiances  just 
below  and  just  above  the  sea  surface,  respectively.  The 
spectral  ratio  given  by  (2)  can  be  expressed  in  terms  of 
/?„(A)  as 

Ty  =  /?.,(A,)//?„(Aj)*£j(A,,  0-i/£j(A^,  0-)  (5) 

or.  for  areas  where  £„(A.  0-)/Z.„(A)  is  only  weakly  depen¬ 
dent  on  wavelength  [Ai«/m.  19791,  as  the  approximation 

r,i  =»  £(A,)/£(Aj)»£^(A,,  0-m/£j(A^.  0-^)  (6) 

since  parameters  associated  with  passing  irradiance  thru  the 
air-sea  interface  are  nonspectral.  For  cloudless  skies.  £^1  A. 
0-L)  can  be  calculated  based  upon  a  number  of  models 
[Justus  and  Paris.  1985;  Bird  and  Riordan.  1986;  Frouin  et 
al..  1989;  Gregg  and  Carder,  1990).  Thus  reflectance  data  of 
either  type  can  be  used  interchangably  in  spectral  ratio  form 
as  long  as  £„(A.  0-)/Z.»(A)  is  relatively  constant  with 
wavelength.  The  implications  of  this  condition  will  be  con¬ 
sidered  in  section  5. 

Our  model  (called  the  “DP  model”)  employs  the  spectral 
ratio  concept  but  uses  a  semianalytical  formulation  of  irra¬ 
diance  reflectance.  R(A),  rather  than  the  empirical  approach 
that  is  manifest  in  (I).  Based  upon  the  work  of  .Morel  and 
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Fig.  1.  Chi  u-specific  absorption  coefficient  at  443  nm  for  pani¬ 
cles  (including  phytoplankton  and  PDP)  versus  (Chi  u|.  Data  are 
from  the  fall  1982  ODEX  cruise,  near-surface  i  ;  <  25  ml  samples 
only. 


Prieiir  [197"’]  and  Gordon  et  al.  [1975],  /?(A)  can  be  e.x- 
pressed  in  terms  of  the  inherent  optical  properties  of  the 
water: 

R(,\)  =0.33(f7,;  -  “al  <’7) 

where  b'  is  the  backscattenng  coefficient,  a  is  the  absorption 
coefficient,  and  the  subscnpts  w.  cdom.  and  p  represent 
water.  CDOM.  and  particles.  The  spectral  dependence  of 
each  term  has  been  left  off  for  convenience.  In  practice. 
<Jcdom  corresponds  to  the  fraction  of  a  given  water  sample 
that  passes  through  a  0.2-^  pore  size  hiter.  while 
corresponds  to  the  fraction  that  is  retained  on  a  Whatman 
GF'F  glass-fiber  filter.  This  approach  ignores  the  absorption 
effects  of  bacteria  and  detntus  in  the  size  range  of  0.2  to 
about  0.7  nm  [Altabet.  1990). 

.A  difficulty  with  the  above  e.xpression  is  that  absorption 
effects  due  to  viable  phytoplankton  pigments,  pheopigmencs. 
detritus,  and  bacteria  are  ail  a  part  of  a,.  The  last  three  of 
these  components,  referred  to  hereafter  as  "particulate 
degradation  products  "  or  PDP.  do  not  contnbute  to  the 
photosynthetic  process,  yet  they  often  play  a  significant  role 
in  determining  the  magn.tude  of  a^,  especially  at  the  shoner 
wavelengths  [Mitchell  and  Kiefer.  I988<i|.  This  effect  is 
apparent  in  Figure  1  where  the  Chi  a-specific  absorption 
coefficient  for  particles  at  443  nm.  0^(443).  is  graphed 
against  [Chi  a]  for  samples  taken  during  the  1982  Optical 
Dynamics  Experiment  (ODEX)  cruise  off  California.  The 
variation  in  a^(443)  is  quite  large,  even  over  a  small  range  of 
[Chi  a]  (0.05  <  [Ch!  a]  <  0.2  mg/m^).  Although  some  of  this 
variation  may  be  due  to  the  package  effect  (see  section  2.4). 
evaluation  of  the  absorption  ratio  a, (435)  ;ap(675)  yielded 
values  as  high  as  7.08  with  a  mean  value  of  4.56.  indicating 
that  significant  detntal  absorption  effects  must  also  be 
present  [see  .Mitchell  and  Kiefer.  1988a:  Roesler  et  al.. 
1989].. 

By  removing  the  PDP  absorption  from  and  regrouping 
it  with  the  CDOM  absorption  term,  an  absorption  coefficient 
"that  represents  all  DP.  including  the  0.2-  to  0.7-/xm  size 
fraction,  is  formed  (ajp  *  flcdom 

which  represents  biologically  utilizable  absorption  ( a >  a^ 
-  apdp).  assuming  that  all  viable  pigments  are  contained  in 


cells  larger  than  the  nominal  GF/F  pore  size  of  0.7  fun.  This 
change  leads  to  the  following  reformulation  of  (7): 

/?(A)  »  0.33(6;,  -  6;)/(rt..  fldp  -  a;[Chl  a])  (8) 

where  a^  is  the  Chi  a-specific  absorption  coefficient  due 
only  to  viable  phytoplankton  pigments.  [Chi  a)  (used  as  a 
surrogate  for  the  concentration  of  all  photosynthetic  pig¬ 
ments)  is  chosen  as  the  index  for  specific  absorption  rather 
than  the  commonly  used  [Chi  a]  [pheo  a]  since  absorption 
due  to  pheopigments  is  included  in  the  DP  absorption  term. 
.As  will  be  shown  in  section  2.3.  this  regrouping  fits  smoothly 
into  the  modeling  scheme  due  to  the  similar  spectral  shapes 
*5f  OjjQni  and  Opjp. 

Tables  are  available  for  both  6'„.  [.Viorel.  1974)  and  a„. 
[Smith  and  Baker.  1981).  Expressions  for  each  remaining 
term  in  (8)  will  be  developed  for  each  of  the  three  wave¬ 
lengths  412.  443.  and  565  nm.  Two  irradiance  ratio  equations 
will  be  formed  from  the  three  irradiance  equations,  forming 
the  basis  for  the  pigment  and  degradation  product  algonthm. 

2.2.  Spectral  Backscattering  Due  to  Particles 

Expressions  for  6),  at  443  nm  and  565  nm  are  adapted  from 
empirical  formulas  developed  by  Cordon  et  al.  [1988].  They 
are 


6"(443)  =  0.0030[Chl  fl]®  " 

(9) 

6;(565)  =  0.0033[Chl  a]® 

(10) 

We  have  substituted  [Chi  a]  for  [C]  and  assumed  that  the 
backscattenng  coefficient  at  565  nm  is  the  same  as  that  at  550 
nm.  the  wavelength  that  is  used  by  Gordon  et  al.  [I988|.  The 
spectral  dependence  of  (9)  and  ( 10)  is  weak  for  the  range  of 
[Chi  u]  values  encountered  in  this  study  ([Chi  a]  s  1.3  mg 
m  .  see  Table  2).  and  particle  backscattenng  will  thus  play 
only  a  minor  role  relative  to  absorption  in  affecting  the 
spectral  irradiance  ratios.  For  waters  where  [Chi  a]  is  of  the 
order  of  10  mg  m or  larger,  this  may  not  be  the  case. 

.An  expression  for  b„  at  412  nm  can  be  developed  using 
arguments  that  parallel  those  of  Cordon  and  Morel  [1983] 
and  Cordon  et  al.  [1988].  We  assume  that  (1)  6),(A)  =  A(A) 
(Chi  where  A  and  B  are  constants;  (2)  at  high 

chlorophyll  I  [Chi  u]  =  20  mg  m"’).  6p(4l2)/6pl443)  =  1.2 
due  to  enhanced  p  gment  absorption,  and  thus  decreased 
backscattering.  at  443  versus  412  nm:  and  (3)  at  low  chloro¬ 
phyll  ({Chi  a]  =  0.05  mg  m"’).  backscattenng  follows  a  A*' 
power  law  [Cordon  and  .Morel.  1983].  yielding  6^(412)/ 
6;,(443)  =  1.0752.  Applying  the  points  in  2  and  3  to  the 
equation  in  I  leads  to 

6^(412)/6;(443)  =  l.l4[Chl  u]®*'*  (II) 

and  substituting  (9)  yields 

6;(412)  *  0.0034[Chl  a]®  -^  (12) 

Different  relationships  are  needed  for  situations  where  back¬ 
scattenng  is  poorly  correlated  with  Chi  a.  within  a  cocco- 
lithophore  bloom  for  example  or  where  suspended  sediments 
or  other  nonalgal  scatterers  are  in  high  proponkm  to  living 
cells. 
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2.3.  Spectral  Absorption  by  Degradation  Products 

Recently,  the  term  CDOM  has  been  employed  to  describe 
that  fraction  of  the  dissolved  organic  matter  which  modifies 
the  color  of  the  water  by  absorption  of  light  ( the  terms  gilvin. 
gelbstoff.  or  yellow  substance  are  often  used  interchangably 
with  CDOM).  The  oceanographic  and  limnological  literature 
is  rife  with  reports  of  the  absorption  effects  of  CDOM 
[James  and  Birge.  193^;  Jerlov,  1955:  Kalle.  1%I.  1966: 
Hojerslev,  1974:  Stuermer.  1975:  Lundgren.  1976:  Kirk. 
1976,  1980:  Bricaud  et  al..  1981:  Carder  et  al.,  19891. 
Similarly.  PDF  often  has  a  significant  effect  on  the  spectral 
absorption  of  ocean  panicles,  especially  at  shon  wave¬ 
lengths  [Kiefer  and  Soohoo.  1982:  Kishino  et  al..  1985: 
\iitcheil  and  Kiefer.  1988a:  Roesler  et  al..  1989).  Both 
CDO.M  and  PDP  are  accounted  for  in  our  degradation 
product  absorption  term. 

In  general,  the  CDOM  found  in  most  of  the  oceans  of  the 
world  can  be  of  either  marine  or  terrigenous  ongin.  How¬ 
ever.  G.  Harvey  (personal  communication,  1991)  has  found 
that  the  molecular  weights  of  the  humic  and  fulvic  acid 
samples  collected  from  the  euphotic  zone  of  oceanic  waters 
in  several  oceans  and  even  in  the  offshore  .Mississippi  River 
plume  (salinity  >  25'ico)  were  less  than  1000.  indicative  of  a 
marine  or  aqueous  rather  than  a  terrigenous  ongin  [Hayase 
and  Tsubota.  1985:  Carder  et  al..  1989).  This  is  important 
because  the  optical  properties  of  manne  humus  are  quite 
different  from  those  of  terrestnal  humus.  In  addition.  Carder 
et  al.  [1989]  argue  that  at  least  for  the  Gulf  of  .Me.xico. 
CDOM  absorption  can  largely  be  accounted  for  by  the 
summed  absorptions  of  manne  humic  acid  (.MHA)  and 
manne  fulvic  acid  (.MF.A).  By  adopting  the  approximation 
‘'cdom  -  ~  degradation  product  absorption  can  be 

expressed  as 

ajp  =  a,,  -  a/*  Updp  (13) 

where  the  subscripts  h  and  /  refer  to  MH.A  and  MF.^. 
respectively.  Note  that  this  approximation  can  only  under¬ 
estimate  CDO.M  absorption. 

Carder  et  al.  [1989|  have  shown  that  spectral  absorption 
curves  for  .MH.A  and  .MF.A  can  be  tit  by  exponential  func¬ 
tions  of  the  form 

a,(A)  =  C,at(450)e^'"‘’°''‘'  (14) 

where  a, (A)  is  the  absorption  coefficient  of  component  .r(.r 
=  h  or/)  at  wavelength  A  in  m"' .  C,  is  the  concentration  of 
.t  in  g  m■^  a*(450)  is  the  specific  absorption  coefficient  .if 
.r  at  450  nm  in  m’  g”' ,  and  S,  is  the  spectral  slope  parameter 
of  x  in  nm''.  They  determined  that  a^(450)  =  0.1304  m’ 
g"',  a^450)  =■  0.0073  m’  g''.  =  0  01 1  nm'',  and  Sf 

=  0.019  nm for  MHA  and  .MF.A  samples  from  the  Gulf  of 
Mexico. 

Furthermore,  it  is  clear  from  the  work  of  Kishino  et  al. 
[1985]  and  Roesler  et  al.  (1989)  that  spectral  absorption 
curves  for  PDP  rather  closely  approximate  the  curve  shapes 
for  marine  CDOM  for  wavelengths  greater  than  412  nm  (sec 
Figure  2).  In  fact,  the  global  mean  spectral  slope  for  deirital 
absorption  data  (Sp^p)  (it  to  an  exponential  curve  like  (14) 
was  found  by  Roesler  et  al.  (1989)  to  be  0.011  nm''.  This 
value  is  the  same  as  the  spectral  slope  parameter  5^  for 
MHA  found  by  Carder  et  al.  (1989),  rendering  the  absorption 
by  the  two  components  spectrally  inseparable.  We  thus 


Fig.  Z.  Relative  absorption  coefficients  versus  wavelength  for 
manne  CDOM  («£*>„  ~  ^  J and  PDP  [from  Kiefer  and 
Soohoo.  19821. 


define  a  weighted  concentration  parameter.  C\.  that  ac¬ 
counts  for  both  MH.A  and  PDP  in  the  same  term: 

C;  =  C*  *  Cpdpa;dp(-»50)/aJ(450)  (15) 

The  prime  indicates  that  the  parameter  is  not  a  true  concen¬ 
tration.  but  one  which  includes  PDP,  which  is  weighted  by 
the  specific  absorption  coefficients. 

U sing  ( 1 5)  and  the  Gulf  of  Mexico  parameters  from  Carder 
et  al.  (1989).  (13)  can  be  rewritten  as 

adp  =  0.1304C;e“  “"'*’®''*’  +  0.0073C/e'’ 

( 16) 

Here  a^p  is  a  function  of  C),  and  C/  for  a  given  wavelength. 
However,  in  order  to  apply  an  algorithm  with  two  irradiance 
ratios,  the  model  is  underspecific  unless  a^p  can  be  param¬ 
eterized  in  terms  of  just  a  single  variable  quantity  (see 
section  2.5).  This  can  be  achieved  by  specifying  a  prion  the 
■MF.A  fraction  of  total  DP.  This  fraction  is  defined  as  f  where 

/  =  QQp  =  C/(C;  -  C/  (17) 

and  C'dp  is  called  the  "weighted  degradation  product  con¬ 
centration."  Again,  the  pnme  is  used  to  indicate  that  the 
weighted  PDP  concentration  is  present  implicitly.  The  full 
equation  for  DP  absorption  is  now 

c‘ip  =  Cip[0 . 1 304(  1  )e 0 1 1 '^50  -  *  ’ 

-  0.0073/?®“'^'^’®'*’]  (18) 

with  C'dp  being  the  sole  vanable  for  a  given  f  and  A.  Note 
that  for  regions  where  absorption  due  to  CDOM  is  high 
relative  to  that  due  to  PDP  (which  is  usually  the  case).  C),  = 
C*  and  C'pp  is  equivalent  to  the  CDOM  concentration. 

It  is  important  to  remember  that  the  parameters  a ^(450). 
aj<450).  5*.  Sf,  and  f.  are  not  universal  and  that  they 
should  be  evaluated  on  a  site-specific  basis.  This  is  particu¬ 
larly  important  in  regions  that  arc  heavily  influenced  by 
terrigenous  and  riverine  CDOM.  since  the  spectral  slopes 
and  especially  the  rnass-specific  absorption  coefficients  for 
humic  and  fulvic  acids  are  dependent  upon  their  molecular 
weights  [Hayase  and  Tsub  yta.  1985:  Carder  et  al..  1989). 
and  since  the  molecular  weights  of  soil-derived  humus  can 
exceed  100.000  (Havase  and  Tsubota.  1985).  On  the  positive 
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side,  XftieUc.  and  Lunye  [1989)  have  demonstrated  that 
certain  very  large  provinces  exist  which  have  stable  or  only 
slowly  changing  bio-optical  parameters. 

2.4.  Spectral  Absorption  hy  Phytoplankton 

In  this  section,  expressions  for  the  Chi  a-specihc  absorp¬ 
tion  coefficient  for  phytoplankton.  n^lA).  are  developed  as 
functions  of  [Chi  a]  for  the  three  wavelengths  of  interest. 
This  task  is  complicated  by  the  complex  bio-optical  charac¬ 
teristics  of  phytoplankton.  Morel  and  Brtcaud[\9i\]  showed 
theoretically  that  n^(A)  is  a  function  of  the  absorption 
coefficient  of  the  cell  material,  the  intracellular  pigment 
concentration,  and  ceil  size.  The  'pigment  package  "  or 
self-shading  effect  that  they  descnbe  results  in  an  inverse 
relation  between  a'^(A)  and  cell  size  for  a  given  internal  cell 
pigment  concentration  and  ceil  matenal  absorption  coeffi¬ 
cient,  and  it  is  manifest  most  significantly  near  the  pigment 
absorption  maximum  of  436  nm.  .Also,  recent  work  has 
shown  that  Ua*-'*  vanes  with  light  [Dubinsky  et  at..  1986: 
.Mitchell  and  Kiefer.  19886]  and  nutnent  adaptation  [Sosik 
and  .Mitchell.  1991  j.  These  variations  are  due  to  the  relative 
internal  concentrations  of  Chi  a  and  photoprotective  pig¬ 
ments.  and  package  effects.  These  lab  studies  suggest  that 
a^IA)  for  surface  waters  may  also  have  a  seasonal  and 
latitudinal  signature,  dependent  on  species  composition  and 
light  and  nutrient  history.  In  order  to  characterize  phyto¬ 
plankton  absorption,  a  conceptual  understanding  is  needed 
of  how  cell  size,  light  history,  and  nutnent  limitation  affect 
the  pigment  packaging. 

To  begin,  consider  two  studies  that  have  pointed  out  the 
correlation  between  cell  size  and  nutnent  availability.  Her- 
bland  et  al.  [19851  found  that  pigments  in  the  submicron 
phytoplankton  size  fraction  dominated  pigment  concentra¬ 
tion  values  for  oligotrophic  surface  waters  in  the  equatonal 
Atlantic,  whereas  they  were  only  a  minor  component  for 
nutrient-  and  chlorophyll-rich  waters  near  the  bottom  of  the 
photic  zone.  Carder  et  al.  [1986]  found  an  analogous  trend  in 
the  onshore-offshore  direction  for  subtropical  waters  off 
Flonda:  that  is.  the  pigments  of  the  nutrient-  and  chloro¬ 
phyll-rich  inshore  waters  were  generally  found  packaged  in 
larger  cells  than  were  the  pigments  of  the  more  oligotrophic 
offshore  waters.  These  findings  suggest  that  for  a  given  light 
domain  (e.g..  season,  latitude).  [Chi  a]  is  roughly  correlated 
with  ensemble  average  cell  size,  which  in  turn  suggests  a 
rough  inverse  correlation  between  [Chi  a)  and  a*^(A). 

.Another  factor  to  consider  is  that  as  [Chi  uj  approaches 
very  low  or  high  values.  u^lA)  approaches  physiological 
limits,  rather  like  asymptotes.  For  decreasing  [Chi  a],  the 
cell  size  and  the  package  effect  decrease  until  a'^tA)  ap¬ 
proaches  the  upper  limit  provided  by  soluble  pigments  which 
is  about  0.10  m"  (mg  Chi  a)*'  at  435  nm  for  high-light, 
high-carotenoid  conditions  [Bricaud  et  al..  1988).  For  in¬ 
creasing  [Chi  a],  a*j(  A)  decreases  because  a  higher  fraction 
of  the  pigments  are  contained  in  larger,  pigment-nch  cells  in 
high-chlorophyll  waters. 

These  two  concepts,  the  inverse  correlation  between 
u^(A)  and  [Chi  a|  and  the  tact  that  there  are  physiological 
upper  and  lower  bounds  on  a ',(  A)  are  used  as  the  basis  for  a 
preliminary  desenption  of  how  a'^lA)  may  behave  as  a 
function  of  [Chi  a].  We  chose  a  hyperbolic  tangent  function 
to  provide  a  curve  shape  that  is  limited  by  asymptotes,  and 
logarithmic  coordinates  are  used  to  deal  with  the  large 
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Fig.  3.  Proposed  functional  relationships  between  Chi 
a-specific  absorption  coefficient  at  443  nm  for  phytoplankton  and 
(Chi  a)  for  subtropical  and  temperate  ocean  regimes.  Plotted  points 
are  ‘detritus-free  '  data  that  are  available  to  corroborate  the  curves. 


dynamic  range  of  [Chi  a]  values  found  in  natural  waten  (see 
Figure  3).  For  the  443-nm  channel  our  expression  has  the 
generic  form 

a;(443)  =  A  exp  {B  tanh  (C  In  ([Chi  a]ID))  (19) 

where  the  parameters  A  and  D  describe  the  ordinate  and 
abscissa  of  the  point  of  symmetry  for  the  curve.  B  describes 
the  asymptotes  of  the  graph,  and  C  desenbes  how  quickly 
the  curve  approaches  the  asymptotes.  Choosing  values  for 
the  parameters  A.  B.  C.  and  D  requires  that  the  curve  be  fit 
to  an  appropnate  set  of  a  *(443)  versus  [Chi  «]  data  points. 
Since  tropical  and  subtropical  phytoplankton  are  typically 
smaller  than  are  cells  from  temperate  and  boreal  waters,  a 
global  empirical  relationship  would  likely  be  less  accurate 
due  to  package  effects  than  would  separate  relationships  for 
different  oceanic  regimes.  Due  to  this  size  variation  and  to 
light  history  considerations,  it  makes  sense  to  define  at  least 
two  separate  regimes,  "subtropical”  and  "temperate.  ’ 

Most  of  the  phytoplankton  absorption  data  used  for  ad¬ 
justing  the  a  *(443)  versus  [Chi  a]  curves  for  both  regimes 
were  determined  by  the  simple  and  inexpensive  GF/F  filter 
pad  transmission  method  described  by  Mitchell  and  Kiefer 
[1984.  19886).  This  method  also  incorporates  absorption  due 
to  particulate  detntus  into  the  measured  absorption  spec¬ 
trum.  so  only  those  data  which  are  reasonably  "detritus- 
free"  have  been  used  (a  methodology  for  separating  the 
absorption  effects  of  phytoplankton  from  those  due  to  detri¬ 
tus  has  been  developed  by  Kishino  et  al.  [1985],  but  unfor¬ 
tunately.  few  data  using  this  methodology  have  been  report¬ 
ed).  Recently.  Stramski  [1990]  reported  that  particulate 
optical  densities  determined  by  this  method  for  two  species 
of  diatoms  changed  rapidly  after  filtration,  presumably  due 
to  degradation  of  chlorophyll  to  pheopigments.  However,  a 
cyanobacterium  and  a  flagellate  species  showed  no  such 
artifacts  when  tested,  and  it  appears  that  absorption  at  443 
nm  decreased  by  only  about  8%  15  min  after  filtration  of  one 
of  the  diatom  cultures.  Since  our  samples  are  measured 
directly  after  filtration,  the  error  introduced  should  be  much 
smaller  than  8%.  especially  if  the  phytoplankton  assemblage 
on  the  filter  pad  is  composed  largely  of  species  that  do  not 
exhibit  the  degradation  effect. 

The  data  sets  for  both  regimes  are  sparse,  but  the  subtrofF- 
ical  data  has  a  greater  range  of  [Chi  u]  values,  and  it  is 
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Fig.  4.  Chi  (i-ipecific  absorption  coelficienl  for  phytoplankton 
versus  wavelength  in  a  high  light-adapted,  oligotrophic  environment 
(25°N.  I55’WI  and  an  eutrophic,  coastal  environment  i  Baja  Califor¬ 
nia).  Picoplanklon  dominate  the  oligotrophic  curve,  where  >65'^  of 
the  Chi  1/  passed  through  a  I.U-Mm  Nuclepore  niter.  High-pressure 
liquid  chromatography  data  for  this  station  vielded  relatively  high 
ratios  of  carotenoid  :  Chi  a  ILjhs  ei  al-.  1990).  This  helps  to  explain 
the  high  peak  at  about  4-kJ  nm  with  values  >0. 10  m*  (mg  Chi  u)  . 
Additional  factors  causing  this  peak  are  the  predominance  of  sub- 
micron  cells,  high  light,  and  low  nutrients,  which  typically  result  in 
low  Chi  a  per  cell.  Detriial  absorption  effects  have  been  removed 
from  the  eutrophic  curve  mathemaiicallv  [Kiefer  and  Soohoo.  19821. 


augmented  on  the  high  [Chi  a]  end  by  data  for  the  Gulf  of 
.Me.xico  that  has  been  derived  from  the  work  of  Carder  et  al. 
(1986).  The  subtropical  data  points  consist  of  oligotrophic 
data  from  the  spring  bloom  north  of  Hawaii  and  eutrophic 
data  from  near  Baja  California  taken  from  Kiefer  and 
Soohoo  (1982).  Laws  et  al.  (1990),  whose  absorption  data  set 
is  the  one  from  which  the  oligotrophic  data  mentioned  above 
derive,  provide  arguments  for  the  detritus-free  nature  of 
chose  data,  while  detrital  effects  have  been  statistically 
removed  from  the  Baja  data.  Figure  3  illustrates  the  pro¬ 
posed  subtropical  curve  shape  (as  well  as  the  proposed 
temperate  curve),  and  it  is  described  by 

ti;(443l  =  0.044 

•  e.xp  { 1 .05  tanh  i  -0.60  in  ([Chi  u]/0.7)))  (20) 

While  there  is  a  need  for  additional  data  points  to  better 
corroborate  the  proposed  curve  shape,  it  does  represent  the 
data  reasonably  well. 

The  high-chlorophyll  data  points  for  the  temperate  curve 
in  Figure  3  were  determined  via  filter  pad  transmission  from 
upwelling  stations  of  the  coastal  transition  zone  cruise  off 
northern  California  in  July  1988.  Only  those  data  which 
exhibited  relatively  low  particulate  absorption  at  400  nm 
relative  to  the  absorption  at  443  nm  are  used.  Detntus-free 
data  for  low-chlorophyll  temperate  waters  have  not  been 
located  however,  precluding  an  attempt  to  empincally  adjust 
the  equation  parameters.  In  light  of  the  limited  data  avail¬ 
able.  we  have  adopted  the  subtropical  curve  shape  and 
adjusted  the  lead  coefficient  downwards  from  0.044  to  0.02  in 
order  to  provide  a  temperate  relationship  for  algorithm 
development  purposes.  Note  that  the  curve  does  appear  to 
provide  a  general  lower  limit  to  the  ODEX  specific  absorp¬ 
tion  data  shown  in  Figure  I.  It  is  this  temperate  regime 
expression  that  is  used  in  the  reflectance  modei. 

The  curves  in  Figure  4  are  used  to  develop  package-effect 


relationships  that  express  a^(4l2)  and  <1^(565)  as  fractions 
of  11^(443).  The  aV  A)  curve  from  a  high-light,  oligotrophic. 
subtropical  environment  where  picoplankton  predominate 
was  obtained  at  26'’N.  liS^W  at  a  depth  of  20  m  on  .April  2. 
1987.  This  curve  was  measured  via  the  filter  pad  transmis¬ 
sion  method  ol  Mitchell  and  Kiefer  (1984.  19886|.  and  it  has 
a  minimaT  package  effect.  The  a“j,(A)  curve  for  high- 
chlorophyll  subtropical  waters  is  the  ensemble-average 
curve  of  Kiefer  and  Soohoo  (1982)  for  coastal  waters  near 
Baja  California  with  the  detrital  absorption  effects  statisti¬ 
cally  removed.  Once  again  the  hyperbolic  tangent  function  in 
logarithmic  coordinates  was  chosen  to  model  these  relation¬ 
ships.  Using  points  from  Figure  4  to  determine  extrema 
yielded  the  following  expressions: 

<i;(4l2) 

=  (0.85  exp  (0.2  tanh  (0.4  in  i[Chl  u]/.6)i}}  u-,(443i 

(21 ) 

u;i565i 

=  {0.2  exp  (0.4  tanh  (0.4  In  ([Chi  a]/.6i)\}  a  Ji443) 

(22) 

U  must  be  emphasized  that  by  no  means  are  the  formula¬ 
tions  developed  in  this  section  meant  to  be  universal.  The 
generic  forms  are  used  in  an  attempt  to  model  the  behavior 
of  a  *(  A)  as  a  function  of  (Chi  a)  based  on  the  theoretical 
concepts  regarding  the  package  effect  discussed  earlier.  In 
light  of  the  limited  data  available,  the  relationships  given  are 
first  approximaiions  that  can  be  tolerated  temporarily  in 
order  to  develop  the  algorithm  as  a  whole.  As  data  sets  of 
adequate  size  emerge,  the  new  relationships  that  result  can 
easily  be  inserted  into  the  model.  Researchers  are  encour¬ 
aged  to  develop  a‘JX)  expressions  for  various  environments 
using  the  improved  qualitative  filter  technique  of  Mitchell 
(1990),  augmented  by  both  detrital  correction  [Kishino  et  al. . 
1985)  and  pheopigment  absorption  correction  [Roesler  et  al. . 
1989). 

2.5.  Pigment  Algorithm 

Inserting  the  appropriate  expressions  for  each  of  the  terms 
in  the  irradiance  reflectance  equations  (equation  (8)  with  A  = 
412.  443.  and  565  nm)  results  in  the  two  irradiance  reflec¬ 
tance  ratios 

/?(412> 

/?(443) 

[f>;(412)  -(•  »;(4l2)][a,.(443)  -  adp(443)  ^  a^(443)] 
[^^(443)  -c  6^(443)][«„.(4I2)  -  ajp(4l2)  a4(4l2)] 

(23) 

/?(443) 

R1565) 

[l>;(443)  •)•  6;(443)][u,,(565)  -  0^^1565)  a^(565l] 

*  [A;(565)  -  <>;(565)][rt,,(443)  +  ajp(443)  +  a^t443)] 

(24) 
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where 

a<ip(A)  =  [0.1304(1  *1 

+  0.0073/’e®®'’‘^’®-*']Cip 
ap(4l2)  =  0.85*’®  -  “''“  ^"0^(443) 

a^^443)  =  o.02e '  '■®  ‘’ '1™ 

ap(565)  =  0.20*-®  •* *  ®  ‘'>d.*(443) 

By  forming  two  irradiance  reflectance  ratios.  /?(412)//?(443) 
and  /?(443)//?(565).  we  have  two  equations  written  in  two 
unknowns.  [Chi  a]  and  CJjp.  Since  analytical  inversion  of 
this  system  of  equations  is  difficult,  nomograms  (Figure  5) 
and  computer  look-up  tables  are  generated  that  relate  paired 
values  of  /?(4l2)//?(443)  and  log(/?(443)//?(565)j  to  varying 
paired  values  of  [Chi  uj  and  Cjp  for  a  given  fiilvic  acid 
fraction  f .  The  right  sides  of  (23)  and  (24)  are  solved  for  a 
two-dimensional  array  of  [Chi  a  \  and  C'jp  values  at  a  fixed/ 
and  the  corresponding  array  of  irradiance  reflectance  ratios 
is  tabulated  in  Table  1.  Forty-six  discrete  concentration 
values  for  each  ocean  color  constituent  were  used  to  con¬ 
struct  the  look-up  tables,  ranging  from  O.Ol  to  3.0  mg  m'^ 
for  [Chi  a]  and  from  0  to  6.0  g  m  for  C'^p.  Three  look-up 
tables  were  constructed,  corresponding  to  fulvic  acid  frac¬ 
tions  of  0.89, 0.92.  and  0.95.  Irradiance  reflectance  ratio  data 
can  be  compared  to  the  calculated  values  in  the  look-up 
tables,  and  estimates  of  [Chi  a]  and  Cjp  can  be  obtained  via 
two-dimensional  linear  interpolation  of  the  input  spectral 
ratios  between  the  tabulated  spectral  ratios.  Quick  estimates 
can  be  made  by  selecting  points  on  the  nomogram  (Figure  5). 

3.  Field  Measurements 

The  optical  field  measurements  used  to  test  this  algorithm 
(referred  to  as  the  "DP  algorithm")  consist  of  irradiance 
reflectance  measurements  taken  in  the  top  attenuation  layer 
( c  <  l/krf(490  nm))  dunng  the  fall  1982  ODEX  cruise  west 
of  California  (sec  Figure  6  for  station  locations).  The  envi¬ 
ronment  included  oligotrophic.  coastal  transition,  and  eu- 
trophic  water  types.  Measurements  of  £„(  A)  and  were 

made  using  the  Biosphencal  Instruments  Bio-Optical  Pro¬ 
filing  System  (BOPS,  described  by  Smith  et  al.  [1984])  for  A 


log  [R(443)  /  R(565)] 

Fig.  3.  Nomogram  of  (Chi  a  \  (mg  m"’)  and  C'jg  (g  m"’)  as  a 
function  of  /J(4i:)/R(443)  and  log(/f(443)//KS65)).  •  C*  - 

C/  Cpip  a^p(450)/a^(450).  Here./  «  0.92. 


Table  I.  summary  of  the  Individual  Scattering  and 
Absorption  Curves 


A 

412 

443 

565 

0.00333 

0.00237 

0.0008.2 

«>p(A) 

0.0034  (Chldl®  -^ 

0.0030  (Chi  ul®“ 

0.0033  (Chi  ^ 

a,(AI 

0.0160 

0.0145 

0.0787 

=  410.  441. 465.  488.  520.  540.  560.  589.  625.  671. 694.  and 
767  nm.  To  minimize  the  effect  of  changing  skies  [Smith  and 
Baker.  1984)  and  ship  shadow  effects  [Gordon.  1985:  Smith 
and  Baker.  1986).  we  only  used  measurements  that  were 
taken  on  sunny  days,  simulating  conditions  appropriate  for 
remote  sensing.  The  BOPS  wavelengths  410.  441.  and  560 
nm  arc  used  as  surrogates  for  the  DP  model  wavelengths  of 
4I2.  443.  and  565  nm,  which  are  the  proposed  wavelengths 
for  Sea-WiFS.  The  appropriate  measured  irradiance  reflec¬ 
tance  ratios  are  given  in  Table  2. 

(Chi  a]  values  are  based  on  in  situ  fluorescence  measured 
at  10  m  depth.  Discrete  chlorophyll  and  pheopigment  deter¬ 
minations  were  made  at  selected  depths  on  each  cast  using  a 
Turner  11  lA  fluorometer.  Millipore  HA  filters  of  pore  size 
0.45  Acm  were  used  to  filter  each  seawater  sample.  These 
filters  dissolve  in  90%  acetone  which  permits  a  no-grind 
chlorophyll  extraction  technique  [Smi7>i  et  al..  1981).  The 
BOPS  provided  measurements  of  continuous  vertical  pig¬ 
ment  fluorescence  which  were  then  calibrated  by  the  dis¬ 
crete  extracted  chlorophyll  values  available  for  each  cast. 

4.  Results 

The  optical  data  set  described  above  has  been  used  for 
input  to  both  the  DP  algorithm  and  to  an  appropriate  Case  1 
algorithm.  First.  R(.K)  is  calculated  from  EJki/Ejik)  for  A 
=  412.  443.  and  565  nm.  For  the  DP  algorithm.  [Chi  a]  and 
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TABLE  2.  Measured  and  Modeled  (Chi  a]  Values  From  OOEX  Cruise.  Fall  1982 


Station 

RtaiO) 

/?t441) 

Chi  a 

(Measured) 

pheo  <1 
(.Measured) 

C 

(Case!) 

Chi  a 
(DP) 

Cjp 

(DP) 

C  jp 

R(441) 

/t(560) 

Chi  a 

9d 

0.922 

1.116 

1.307 

0.42 

1.400 

1.909 

3.355 

9u 

0.903 

1.129 

1.311 

0.42 

1.372 

1.337 

3.846 

2.93 

21d 

0.%5 

2.877 

■0.130 

0.06 

0  250 

0.191 

1.419 

10.91 

21u 

0.926 

2.839 

0.121 

0.05 

0.256 

0.109 

1.831 

15.13 

:3d 

0.988 

3.973 

0.128 

0.04 

0.139 

0.104 

0.959 

'.49 

23  u 

1.007 

4.247 

0.128 

0.04 

0.123 

0.106 

0.808 

6.31 

24d 

0.987 

4.263 

0.112 

0.03 

0.122 

0.084 

0.911 

8.13 

29d 

0.917 

2.156 

0.242 

0.10 

0.423 

0.201 

2.390 

9.87 

30d 

0.992 

3.492 

0.136 

0.05 

0.176 

0.152 

1.046 

'  69 

30u 

0.998 

3.520 

0.139 

0.05 

0.173 

0.161 

1.000 

'.19 

36d 

1.020 

4.295 

0.122 

0.04 

0.121 

0.116 

0.'41 

6.07 

47d 

1.102 

6.029 

0.085 

0.01 

0.065 

0.088 

0.303 

''.56 

77.  Id 

1.109 

6.522 

0.074 

0.03 

0.056 

0.074 

0.259 

3.50 

77.:d 

1.173 

6.659 

0.074 

0.03 

0.054 

0.099 

0.130 

1.'6 

84u 

1.142 

6.328 

0.102 

0.06 

0.060 

0  097 

0.200 

1.96 

93u 

1.101 

6.087 

0.091 

0.09 

0.064 

0.085 

0.303 

3.33 

94d 

1.141 

5.999 

0.105 

0.08 

0.066 

O.lll 

0.221 

2.10 

95d 

1.097 

6.135 

0.113 

0.02 

0.063 

0.081 

0.309 

2. '3 

I04u 

1.131 

6.142 

0.108 

0.01 

0.063 

0.099 

0.234 

2.16 

170. 2d 

1.045 

4.317 

0.140 

0.04 

0.119 

0.141 

0.637 

4.54 

I70.2U 

1.046 

4.374 

0.142 

0.04 

0.117 

0.138 

0.624 

4.39 

I74d 

0.937 

3.327 

0.158 

0.04 

0.192 

0.084 

1.508 

9  54 

|75.1d 

0.946 

2.963 

0  253 

0.00 

0.237 

0.135 

1.564 

6.18 

176.  Id 

0.901 

1.825 

0.308 

0.1 1 

0.572 

0.238 

3.010 

9.77 

176.  lu 

0.897 

1.860 

0.299 

0.11 

0.553 

0.202 

3.084 

10.31 

176.2d 

0.900 

1.589 

0.328 

0.12 

0.736 

0.360 

3.329 

10.15 

From  lett  to  right  the  columns  represent  the  OOEX  station  number,  measured  irradiance  reflectance  ratios,  measured  [Chi  a]  and  (pheo 
u].  (Cl  calculated  with  the  case  I  algorithm.  (Chi  u|  and  C^p  calculated  with  the  DP  algotithm.  and  the  ratio  of  (modeled  Cjp):  (measured 
Chi  a).  [Chi  j|  and  [pheo  a\  are  in  mg  m  and  Cjp  is  approximately  equal  to  CDO.M  concentration  in  g  m'V 


C'jp  are  calculated  from  R{\)  values  as  described  in  section 
2.5.  and  for  the  case  I  algonthm.  [C]  is  calculated  by 

[C]=  I.71[/?(440)//?(560i]‘' (25) 

which  is  described  in  further  detail  by  Gordon  and  \forel 
(1983).  (Chi  a|  estimated  from  the  DP  algorithm  and  (C) 
derived  from  the  Case  1  algorithm  are  compared  with 
measured  (Chi  a\.  and  the  performance  of  the  two  algorithms 
are  tested.  Only  the  accuracy  of  the  Chi  a  ponion  of  the  DP 
algorithm  could  be  tested  since  measurements  of  DP  absorp¬ 
tion  are  unavailable  for  comparison  with  the  model  estimates 
of  Cjp.  The  results  are  given  in  Table  2, 

Before  employing  the  DP  algorithm  it  is  necessary  to 
specify/’,  the  fulvic  acid  fraction  of  total  CDOM.  The  DP 
algonthm  was  run  three  times  with/  held  constant  for  each 
run  at  values  of  0.89,  0.92.  and  0.95.  The  total  average  error 
in  calculated  (Chi  a)  versus  measured  (Chi  a\  was  the  lowest 
for/’  =  0.92.  and  it  is  assumed  that  0.92  is  the  '’regional 
average  '/’.  It  is  the  results  from  this  run  that  are  shown  in 
Table  2.  Used  in  this  manner,  the  DP  algorithm  cannot  be 
purely  predictive  without  some  prior  knowledge  of/.  For¬ 
tunately.  as  will  be  seen  in  section  5.  this  problem  should  not 
be  large. 

The  results  are  depicted  graphically  in  Figure  7.  Figure  7a 
shows  (Chi  a\  derived  from  the  DP  algorithm  and  (Cl 
derived  from  the  case  I  algorithm  versus  measured  (Chi  a] 
for  all  26  stations.  High  Chi  a  points  ((Chi  «|  >  I)  are  well 
described  by  the  case  I  algorithm,  whereas  the  DP  algorithm 
estimates  are  less  accurate.  Two  observations  may  help  to 
explain  this.  First.  Table  2  reveals  that  the  ratio  (modeled 
Cjp) ;  (measured  (Chi  aj)  is  relatively  low  for  these  waters. 


indicating  that  pigment  absorption  is  the  dominant  influence 
on  the  f?(443) ;  /;i565)  reflectance  ratio:  the  case  I  algorithm 
works  well  when  noncovarying  optical  components  (e.g.. 
CDOM)  are  not  high.  Second,  the  nomogram  (Figure  5) 
shows  that  the  DP  algorithm  becomes  increasingly  sensitive 
to  changes  in  the  input  reflectance  ratios  at  higher  Chi  a  and 
C|jp.  Thus  the  DP  algorithm  will  be  less  accurate  at  high  Chi 
a  than  in  regions  of  medium-to-low  Chi  a  ((Chi  a]  <  I). 

Figure  7b  is  the  same  as  Figure  7a  except  (hat  only  the 
medium-tolow  Chi  a  points  are  plotted.  At  this  greater 
resolution  the  case  I  algorithm  can  be  characterized  by  two 
distinct  regions:  a  low  Chi  a  region  in  which  (Chi  a\ 
estimates  are  consistently  lower  than  the  measured  values 
and  a  medium  Chi  a  region  in  which  they  are  consistently 
higher  than  the  measured  values.  The  DP  algonthm  predicts 
(Chi  a]  more  accurately  in  these  regions,  showing  little  or  no 
net  bias  at  low  chlorophyll  and  only  a  slight  low  bias  in  the 
medium  region. 

Another  way  to  present  these  data  is  to  plot  the  percent 
error  in  (Chi  a|  estimates  for  the  two  methods  versus  the 
C’jp  :(Chl  a]  ratio  (see  Figure  8).  The  mean  fractional  error 
for  all  stations  was  i38%  for  the  case  I  algorithm  and  r  18% 
for  the  DP  model.  The  highest  errors,  up  to  13396.  occur  in 
using  the  case  I  model  for  waters  with  the  highest  C’pp  :  (Chi 
aj. 

The  percent  error  for  the  case  I  algorithm  in  Figure  8 
correlates  well  with  the  C^p  ;(Chl  a]  ratio,  suggesting  that 
this  ratio  can  be  used  to  indicate  which  waters  are  well 
described  by  the  case  I  algorithm  and  which  waters  are  not. 
For  this  particular  cruise  at  least,  an  operational  means  of 
establishing  the  boundaries  between  true  case  I  and  DP-rich 
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waters  is  to  designate  the  case  1  waters  as  those  for  which 
C'(jp  :[Chl  I/I  <  7  g  mg"'.  Waters  with  high  Cjp  values  are 
not  necessarily  DP-rich  by  this  definition,  notably  stations  Od 
and  9u.  Subdividing  the  study  area  in  such  a  way  is  consis¬ 
tent  with  the  concept  of  “bio-optical  provinces"  \PUtit  and 
Sathyendranath,  1988:  Mueller  and  Laniie.  19891.  which 
states  that  the  absorption  and  scattering  properties  of  sea¬ 
water  vary  from  one  geographical  region  to  the  next  and  with 
respect  to  time  (e.g..  seasonal  variation). 

Referring  to  the  station  location  map  ( Figure  6)  and  to  the 
tabulated  Cjp  :  [Chi  a]  values  (Table  2).  it  can  be  seen  that  in 
general,  the  DP-rich  waters  are  found  seaward  of  station  9 
and  landward  of  station  36.  Thus  both  the  pigment-nch 
coastal  upwelling  waters  (represented  here  solely  by  station 
9)  and  those  waters  found  beyond  the  coastal  transition 
waters  more  than  400  km  olfshore  are  case  1  waters  using 
this  criterion,  while  those  found  in  the  transition  waters 
between  the  coastal  and  oflFshore  waters  are  DP-rich.  The 
only  exceptions  to  this  geographic  scheme  are  stations  23u 
and  175. Id.  which  are  within  the  oceanographic  DP-rich 
boundaries  but  which  barely  miss  the  DP-rich  catena.  Given 
the  rather  chaotic  nature  of  currents  oflF  the  west  coast  of 
California  [see  Abbott  et  al..  1990],  this  is  well  within  the 
scope  of  expected  patchiness.  Also,  the  application  of  this 
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Fig.  7.  Plot  of  [Chi  u|  derived  from  the  DP  model  and  |C| 
derived  from  the  case  I  model  ([C]  »  1.71  [/?(440i//f(}60)r'  *-i 
versus  measured  [Chi  (i|.  Figure  7a  includes  all  data  points,  while 
Figure  7b  is  a  magnification  of  the  lower  conceniraiior  range.  The 
dashed  line  is  the  1:1  slope  line. 
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Fig.  8.  Percent  error  I  = 'modeled  [Chi  ((1  measured  iChl  «!-  1 1 
X  l(X))  for  the  DP  model  and  the  case  I  model  as  a  function  of  ihe 
ratio  of  modeled  Cjp  to  measured  [Chi  u].  The  long  dashed  line  is 
the  zero  line.  The  shon  dashed  lines  are  linear  regression  lines  for 
the  two  data  sets:  )'  =  -  1 .53“.V  -  a  88lr-  =  0.05.  n  =  261  for  the 
DP  model  data,  and  Y  =  I2.8“.V  -  64,7(r-  =  0.80.  n  =  26)  for  the 
case  I  model  data. 


criterion  depends  on  the  estimated  value  of  Cjp.  introducing 
another  modicum  of  imprecision. 

For  the  1 1  DP-rich  stations,  the  mean  error  using  the  case 
I  algorithm  for  [Chi  al  was  =61%.  while  that  for  the  DP 
algonthm  was  =23%.  The  net  positive  bias  in  the  case  I 
estimates  for  these  waters  will  yield  a  corresponding  infla¬ 
tion  of  any  primary  productivity  estimates  that  are  based  on 
such  values,  while  the  DP  algorithm  provided  only  a  slight 
negative  bias  for  these  same  waters.  For  the  case  I  waters 
(15  stations)  the  case  1  algorithm  had  a  mean  fractional  error 
of  =22%  with  a  net  negative  bias  and  the  DP  algorithm  had 
=  14%  error  with  little  or  no  net  bias.  The  error  percentages 
are  summarized  in  Table  3. 

5.  Discussion 

5.1.  Comparison  of  the  DP  Algorithm 
Versus  the  Case  I  Algorithm 

The  most  significant  result  of  this  study  is  the  considerable 
improvement  in  [Chi  «]  retrieval  accuracy  of  the  DP  algo¬ 
rithm  over  the  case  1  algorithm  1 23%  error  versus  6 1  %  error) 
for  those  waters  which  our  model  estimates  to  be  DP-rich 
and  where  [Chi  a)  s  1.3  mg  m*’.  .Although  the  case  1 
algorithm  estimates  [C].  the  sum  of  [Chi  //).  and  [pheo  </|. 
Table  2  demonstrates  that  [C]  and  [Chi  n]  differ  by  only 
10-40%  and  the  presence  of  [pheo  </]  can  only  explain  a 
minor  fraction  of  the  case  I  [Chi  </)  error  for  the  DP-rich 


TABLE  3.  .Mean  Fractional  Error  in  [Chi  o|  Estimates  for  Both 
Case  I  and  DP  .Algorithms 


Data  Subset 

.Algorithm.  ^ 

Case  1 

DP 

n 

.All  Stations 

38 

18 

26 

Cjpi  (Chi  «l  <  7  g  mg  ■' 

14 

1; 

Cjp:  jChl  (ij  >  7  g  mg  ■■ 

61 

23 

II 

Percent  error  is  [I  modeled  [Chi  al)((  measured  (Chi  u|)  -  l|  *  100. 
»  IS  number  of  stations. 


;0.608 
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sTations.  Some  other  color  constituent  must  be  causing  the 
large  overestimates  of  [Chi  a|  generated  by  the  case  1 
algorithm  presented  here.  Since  there  is  no  significant  source 
of  rivenne  or  terrestrial  CDOM  in  our  study  area  and  since 
detrital  particles  are  likely  to  settle  out  quickly  las  Abbott  et 
al.  [1990]  found  off  northern  California),  marine  CDOM  is 
the  only  reasonable  component  of  the  DP  ensemble  that 
remains  to  explain  the  optical  effects  that  we  observe.  This 
interpretation  contradicts  the  assertions  that  marine  CDOM 
always  covaries  with  algal  pigments  or  that  it's  measurable 
influence  is  everywhere  weak  [Morel  and  Prienr.  1977; 
Bricaud  et  al..  1981;  Gordon  and  Morel.  1983).  These 
assertions  are  further  undermined  by  two  previously  men¬ 
tioned  studies  in  which  absorption  due  to  CDOM  increased 
with  time  in  the  same  water  mass  {Peacock  et  al.  [1988)  and 
Carder  et  al.  [1989);  see  section  I). 

In  addition.  [Chi  a]  estimates  from  the  DP  algonthm  are 
better  suited  for  use  with  primary  productivity  models  than 
are  those  of  the  case  1  algorithm.  First.  [Chi  u|  estimates 
from  the  DP  algonthm  are  relatively  free  of  bias  with  respect 
to  the  measured  values,  whereas  [C)  estimates  from  the  case 
1  algorithm  have  considerable  high  bias  in  the  DP-rich  areas 
of  the  study  region  and  low  bias  in  the  case  1  areas  (see 
Figure  8).  These  biases  will  lead  to  systematic  regional 
errors  in  primary  productivity  estimates.  Second,  pnmary 
productivity  models  that  are  based  on  physiological  consid¬ 
erations  need  to  be  driven  by  photon  absorption  by  viable 
pigments.  The  DP  algorithm  explicitly  estimates  [Chi  a|.  the 
major  photosynthetic  pigment,  whereas  the  pigment  esti¬ 
mates  of  the  various  case  1  algorithms  (including  Gordon 
and  .Morel  [1983),  Gordon  et  al.  (1988).  and  Morel  [1988]) 
inherently  include  pheopigments  and  other  degradation 
products  which  do  not  contribute  to  the  photosynthetic 
process. 

5.2.  Regression  Analyses 

The  regression  analyses  performed  on  the  data  depicted  in 
Figure  8  indicate  that  80%  of  the  variance  in  error  values  for 
the  case  I  [Chi  a]  algorithm  could  be  accounted  for  by 
changes  in  the  Cjp:[Chl,al  ratio,  while  only  5%  of  the 
variance  in  error  values  for  the  DP  [Chi  a\  algonthm  error 
could  be  so  explained.  Caution  must  be  used  in  interpreting 
these  results,  not  only  because  the  Cip  values  are  estimates 
and  not  actual  measurements  but  also  because  their  values 
are  inherently  linked  to  the  DP  algorithm  [Chi  a]  values  with 
which  they  are  being  compared.  Prieur  and  Sathyendranath 
[1981]  used  ridge  regression  to  overcome  this  problem  of 
correlation  between  independent  variables.  Since  our  only 
concern  here  is  that  the  covariance  between  the  case  I  error 
and  Cjp  .[Chl  a]  appears  to  be  high,  suggesting  that  a 
non-case  I.  site-specific  [Chi  a\  algorithm  of  the  type  de¬ 
scribed  by  equation  (1)  (i.e.,  a  CZCS-type  algorithm)  can  be 
developed,  the  correlation  problem  was  not  investigated 
further.  .A  linear  regression  of  log  (measured  [Chi  a])  versus 
log(/?(440)/f?(560))  {r-  =  0.84.  n  =  26)  results  in  the 
following  "case  I  case  2"  algorithm: 

[Chi  a]  =  0.80[/?(440)//?(560l]  ■'  (26) 

in  comparison.  .Morel  [1980]  describes  a  similar  'case  I 
case  2"  reflectance  ratio  algorithm,  given  by 


(27) 


(r*  =  0.76.  n  »  90;  see  Gordon  and  Morel  [1983]  for  more 
details),  which  has  a  similar  slope  but  a  lead  coefficient 
which  is  twice  as  large  as  the  algorithm  for  the  data  in  this 
paper.  This  difference  can  be  accounted  for  in  part  by  the 
fact  that  Morel’s  algorithm  is  for  [C].  which  contains 
pheopigments.  whereas  the  ODEX  data  do  not.  and  that 
Morel's  data  set  includes  largely  tropical  and  subtropical 
data  where  there  are  likely  larger  Chi  o -specific  absorption 
coefficients  at  440  nm  than  there  are  for  the  more  temperate 
ODEX  data  (see  Figure  3).  The  algorithm  expressed  by  (26) 
provides  a  means  to  estimate  the  viable  pigments  for  the 
low-light,  late  fall  time  penod  in  the  ODEX  region  using  a 
CZCS-type  algonthm. 


5.3.  Variability  off 

The  results  of  the  DP  algorithm  assume  that  f  =  0.92. 
and  some  of  the  error  in  our  DP  model  results  may  be 
explainable  in  terms  of  a  variable/.  However,  iterative 
trials  of  the  DP  algorithm  in  which/  was  allowed  to  vary  by 
-0.03  about  the  "regional  average  "  of  0.92  for  each  station 
yielded  little  or  no  improvement  in  the  mean  fractional  error 
for  [Chi  n]  estimates,  so  using  a  single  value  of  /  for  the 
entire  study  region  is  probably  not  a  bad  approximation. 

However,  on  a  larger  spatial  scale./  varies  significantly. 
Carder  et  al.  [1989]  evaluated  the  spectral  slope  values  of 
Bricaud  et  al.  [1981]  in  terms  of  MHA  and  MFA  absorption 
and  found  that  their  global  average  value  (5  =  0.014  nm~' 
from  370  nm  to  440  nm)  was  consistent  with  an  /  value  of 
0.89.  Similarly,  the  spectral  slope  for  the  western  north 
.Atlantic  data  of  Topllss  et  al.  [1989]  (S  *  0.013  nm"')  is 
consistent  with  an  t'  value  of  0.84.  Analyzing  1 1  published 
values  of  /  taken  from  the  Gulf  of  Mexico  [Harvey  et  al.. 
1983;  Carder  et  at..  1989]  yields  a  mean  of  0.878  and  a 
sample  standard  deviation  of  0.098.  This  e.xercise  indicates 
that  selection  of  an  appropriate  mean  value  of  /  for  use  in 
the  DP  model  is  site-dependent.  For  river  plumes  and  areas 
heavily  influenced  by  terrestrial  runoff,  the  entire  a^  term 
needs  reassessment,  since  the  geochemistry  of  CDOM  is 
different  for  terrestrial  and  riverine  versus  marine  ecosys¬ 
tems.  resulting  in  significant  variation  in  the  humic/fulvic 
ratio,  the  mass-specific  absorption  coefficients,  and  the 
CDOM  molecular  weights. 


5.4.  DOC  Estimate 

The  concentration  parameter  C’^p  can  provide  a  rough 
measure  of  dissolved  organic  carbon  (DOC)  if  the  factors  in 
the  following  equation  can  be  estimated: 


Cip  (C*  +  Cf  DOC„v 


(28) 


where  DOCy,  is  the  total  DOC  measurable  by  the  ultraviolet 
oxidation  technique  [Armstrong  et  al..  1966]  and  DOCp,  is 
the  total  DOC  measurable  by  the  new  platinum-catalyst 
method  [Sukiiki  et  al..  1985;  Sugimitra  and  Suzuki.  1988]. 
The  factor  (C*  r  C/)/C'jp  corrects  C’^p  for  PDF  and  can  be 
approximated  by  the  easily  measurable  ratio  fl.A--(4IO)/ 
adpl^lO).  The  factor  DOCu»/(C),  Cf  adjusts  for  the  fact 
that  DOC  does  not  consist  only  of  humic  and  fiilvic  acids, 
and  the  final  factor  provides  a  correction  for  the  platinum- 
catalyst  method  versus  the  older  ultraviolet  technique.  .A 
DOCu»/(C),  *  Cf)  value  of  about  3.5  is  provided  by  Har\'ey 


[C]  -  l.62[f?(440)/i?(560)] 
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et  al.  (19831  for  both  oligotrophic  waters  ( ;  =  20  m).  anil  for 
a  low-chlorophyll.  high-CDOM  (i.e.,  DP-rich)  station,  both 
in  the  Gulf  of  Mexico.  Dmjfel  et  al.  [1989]  found  that  the 
p'latinum-catalyst  method  and  the  (JV  oxidation  technique 
yielded  DOC  values  of  about  210  and  90  mg  at  DOC  m*'. 
respectively,  in  the  top  100  m  at  Bl’N.  159°W.  which 
provides  a  DOCp,/DOCuv  ratio  of  2.3.  If  (C^  *  C/)/Cjp  is 
near  1 .0.  these  factors  yield  a  mean  estimated  DOC  concen¬ 
tration  of  1.93  g  DOC  m'’  for  the  cluster  of  seven  olig- 
otrophic  stations  in  this  study  that  are  west  of  140°W  and 
near  33°-34°N.  Using  these  same  factors,  the  12  DP-nch 
stations  in  this  study  have  a  mean  estimated  DOC  concen¬ 
tration  of  14.8  g  DOC  m  The  oligotrophic  ODEX  data  and 
the  Druffel  et  al.  [1989]  data  mentioned  above  are  both  taken 
from  the  North  Pacific  gyre,  so  the  proximity  of  our  estimate 
(1.93  g  DOC  m"^)  to  those  of  Druffel  et  al.  (I989|  (2.52  g 
DOC  m'^)  is  reassuring.  However,  the  steps  between  the 
modeled  C'jp  and  derived  DOC  estimates  contain  general 
approximations  that  have  not  been  tested  at  the  study  site. 
Care  in  interpreting  Cjp  in  terms  of  DOC  should  be  e.xer- 
cised  at  all  study  sites,  especially  those  in  coastal  waters 
where  DOCuv  lC,  —  C^)  ratios  can  be  quite  variable  (see 
Harvey  et  al.,  1983]. 

5.5.  Optical  Effects  on  Water-Leaving  Radiance 

The  DP  model  presented  here  is  based  on  irradiance 
reflectance  and  is  thus  well-suited  for  use  with  moored  or 
free-floating  buoys  with  instrumentation  located  at  5-10  m 
depth  (depths  typical  of  the  data  set  utilized  in  this  study). 
However,  in  judging  the  applicability  of  this  model  for  use 
with  water-leaving  radiance  data,  additional  concerns  must 
be  addressed.  In  making  the  transition  from  upwelling  irra¬ 
diance.  £„(A.  10  m).  to  water-leaving  radiance.  £,„(A). 
consideration  must  be  given  to  changes  in  /?(A)  with  depth, 
spectral  variations  in  £„(A.  0-)/£„(A.  0-),  and  transpec- 
tral  effects  due  to  water  Raman  scattering  and  solar  stimu¬ 
lated  fluorescence  from  CDOM. 

Monte  Carlo  simulations  by  Kirk  (1984)  show  that  Rik) 
changes  very  little  with  depth  for  inelastic  scattering,  so  this 
effect  should  be  negligible.  However,  although  the  wave¬ 
length  dependence  of  the  radiance  distnbution  parameter 
£„(A.  0-)/£„(A.  0-)  is  typically  taken  to  be  spectrally 
quite  flat,  it  apparently  becomes  more  variable  for  A  <  440 
nm  and  A  >  550  nm  (Ausrifi.  1979;  R.  Smith,  unpublished 
Sargasso  Sea  data.  1991;  C.  Davis,  unpublished  California 
Current  data.  1991).  This  problem  cannot  be  so  easily 
ignored. 

In  addition,  recent  measurements  and  model  results  indi¬ 
cate  that  for  oligotrophic  environments,  Raman  inelastic 
scattering  is  not  negligible  [Marshall  and  Smith.  1990;  Stavn. 
1990:  Peacock  et  a!..  1990!  and  that  for  CDOM-rich  environ¬ 
ments.  solar  stimulated  fluorescence  of  the  CDOM  may 
likewise  not  be  negligible  [Hawes  and  Carder.  1990;  Pea¬ 
cock  et  al..  1990).  Contributions  to  the  photon  flux  by  these 
sourcelike  phenomena  will  exert  a  maximum  effect  in  the 
blue-green  wavelengths  near  the  sea  surface  where  UV  and 
near-UV  light  is  available  to  stimulate  emission  of  blue  and 
green  photons.  This  is  so  because  both  the  Raman  scattering 
cross  section  and  CDOM  absorption  increase  with  decreas¬ 
ing  wavelength  and  because  the  availability  of  UV  light  for 
excitation  decreases  rapidly  with  depth. 

An  understanding  of  the  effects  of  these  trnnspectral 


phenomena  on  spectral  ratios  of  £k^(A)  is  necessary  in  order 
to  modify  the  model  and  algorithm  appropriately  for  use  with 
satellite-derived  £„(A)  data.  Peacock  et  al.  (1990)  estimated 
these  effects  to  be  of  second-order  importance,  affecting 
/-w(A)  by  less  than  15%.  If  this  is  indeed  the  case,  the 
modifications  necessary  should  not  be  major. 

5.6.  Conclusion 

For  the  California  Current  region  where  [Chi  a  ]  s  1 . 3  mg 
m'\  the  confounding  effects  of  absorption  due  to  primary 
productivity  degradation  products  on  remote  assessment  of 
chlorophyll  a  can  be  removed  to  a  large  extent  by  the 
addition  of  a  blue  spectral  channel  at  about  412  nm  to  the 
nominal  suite  of  CZCS-like  channels  expected  on  Sea-WiFS 
and  other  future  ocean  color  scanners.  It  is  anticipated  that 
DP  models  such  as  the  one  described  here  will  be  somewhat 
site-  and  perhaps  season-specific  and  that  they  will  signifi¬ 
cantly  enhance  the  accuracy  of  remotely  sensed  estimates  of 
(Chi  a]  for  waters  rich  in  marine  productivity  degradation 
products.  .A  noteworthy  feature  of  the  DP-model  is  that  the 
individual  absorption  and  scattering  parameters  can  be  cus¬ 
tomized  for  a  given  province.  A  particularly  important  task 
is  to  develop  empirical  descriptions  of  the  vanation  in  a^(  A) 
with  respect  to  (Chi  a),  which  is  due  primarily  to  variations 
in  phytoplankton  species  composition  and  light-  and  nutri¬ 
ent-histories.  Presently,  the  technique  best  suited  to  deter¬ 
mine  these  relationships  is  the  improved  quantitative  filter 
technique  of  Mitchell  [\990].  augmented  by  both  the  detrital 
correction  method  of  Kishino  et  al.  (1985)  and  the  pheopig- 
ment  absorption  correction  of  Roesler  et  al.  (1989). 

Notation 

A  wavelength  of  light,  nm. 
fljIA)  absorption  coefficient  for  any  component 
.T.  at  wavelength  A.  m"'  (subscripts  <6. 
cdom.  dp./,  h.  p.  pdp,  and  >v  refer  to 
phytoplankton.  CDOM.  DP.  MF.A.  MHA. 
particles.  PDP.  and  water). 
a*(A)  mass-specific  absorption  coefficient  for 
any  component  .r  at  wavelength  A.  m’ 
g"'  (subscripts  described  above). 
iJ^lA)  Chi  a-specific  absorption  coefficients  for 
phytoplankton  and  particles  at  wavelength 
A.  m'  [mg  Chi  a]  . 

b'^(k).  6',  (A)  backscattering  coefficients  for  particles 
and  water  at  wavelength  A.  m*'. 

(Cl  concentration  of  chlorophyll  a  - 
pheophytin  a.  mg  m‘^ 

C,  concentration  of  component  .r.  g  m  ’’ 
(subscripts  described  above). 

^dp-Ci,  weighted  concentration  parameter  for  DP 
and  MHA.  g  m'^  (the  prime  indicates 
that  the  values  are  weighted  by  the 
presence  of  PDP  via  the  specific 
absorption  coefficients  (see  equation  (15)). 
DOCp,.  DOCuy  dissolved  organic  carbon,  as  determined 
by  the  platinum  catalyst  method  and  by 
the  ultraviolet  oxidation  method,  g  m~^ 
£rf(A.  :).  £„(A.z)  downwelling  and  upwelling 

irradiances  at  wavelength  A  and  depth 
W  m‘‘  nm"'. 


f 


f  fulvic  acid  fraction  of  DP  (the  prime 
indicates  that  the  fraction  is  weighted  by 
the  presence  of  PDF  via  the  specific 
absorption  coefficients  (see  equations  1 15) 
and  (17)). 

A;^i490)  diffuse  attenuation  coefficient  for 

down  welling  irradiance  at  490  nm.  m  . 
L„(K)  subsurface  upwelling  radiance  at 
wavelength  A.  W  m"-  sr*'  nm''. 

Z.„.(A)  water-leaving  radiance  at  wavelength  A.  W 
m*’  sr"'  nm*'. 
n  number  of  samples. 
r,j  ratio  of  water-leaving  radiances  at 
wavelengths  A,  and  A,. 
r-  correlation  coefficient  of  regression. 

/?(A»  irradiance  refiectance  at  wavelength  A. 
/?rs(Ai  remote  sensing  refiectance  at  wavelength 
A. 

5,  spectral  slopes  for  the  exponential 

expressions  that  describe  absorption  due 
any  component  x.  nm*'  (subscripts /,  h. 
and  pdp  refer  to  MF.A.  .VIHA.  and  PDP). 
c  depth,  m. 
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